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ABSTRACT 

Context. Low-metallicity star-forming dwarf galaxies are prime targets to understand the chemical enrichment of the interstellar 
medium (ISM). The H i region contains the bulk of the mass in BCDs, and it provides important constraints on the dispersal and 
mixing of heavy elements released by successive star-formation episodes. The metallicity of the H i region is also a critical parameter 
to investigate the future star-formation history, as metals provide most of the gas cooling that will facilitate and sustain star formation. 
Aims. Our primary objective is to study the enrichment of the H i region and the interplay between star-formation history and metal- 
Ucity evolution. Our secondary objective is to constrain the spatial- and time-scales over which the H i and H ii regions are enriched, 
and the mass range of stars responsible for the heavy element production. Finally, we aim to examine the gas heating and cooling 
mechanisms in the H i region. 

Methods. We observed the most metal-poor star-forming galaxy in the Local Universe, I Zw 18, with the Cosmic Origin Spectrograph 
onboard Hubble. The abundances in the neutral gas are derived from far-ultraviolet absorption-lines (H i, C ii, C ii*, N i, O i, ...) 
and are compared to the abundances in the H ii region. Models are constructed to calculate the ionization structure and the thermal 
processes. We investigate the gas cooling in the H i region through physical diagnostics drawn from the fine-structure level of C"^. 
Results. We find that H i region abundances are lower by a factor of ~ 2 as compared to the H ii region. There is no differential 
depletion on dust between the H i and H n region. Using sulfur as a metallicity tracer, we calculate a metallicity of 1/46 Zq (vs. 
1/31 Z© in the H ii region). From the study of the C/O, [O/Fe], and N/O abundance ratios, we propose that C, N, O, and Fe are mainly 
produced in massive stars. We argue that the H i envelope may contain pockets of pristine gas with a metallicity essentially null. 
Finally, we derive the physical conditions in the H i region by investigating the C n* absorption Une. The cooling rate derived from 
C n* is consistent with collisions with H" atoms in the diffuse neutral gas. We calculate the star-formation rate from the C n* cooling 
rate assuming that photoelectric effect on dust is the dominant gas heating mechanism. Our determination is in good agreement with 
the values in the literature if we assume a low dust-to-gas ratio (~ 2000 times lower than the Milky Way value). 
Conclusions. 

Key words. Galaxies: abundances, HII regions. Galaxies: individual: IZwl8, Galaxies: dwarf. Galaxies: ISM, Galaxies: star forma- 
tion. Galaxies: evolution 



1. Introduction 

Numerical simulations by Cen & Ostriker (1999), predicting 
the evolution of the metal content of the Universe, show that 
metallicity is a stronger function of density than age. It is there- 
fore plausible that gas-rich dwarf galaxies in the Local Universe 
may remain in an early stage of chemical enrichment, possi- 
bly preserving pockets of pristine gas (Kunth & Sargent 1986; 
Kunth et al. 1994). Blue Compact Dwarf (BCD) galaxies rep- 
resent a major testbed as they include the most-metal poor 
star-forming galaxies known (Kunth & Ostlin 2000). Among 
these, IZwl8 (Mrkll6) continues to attract considerable in- 
terest and to feed intense debates in extragalactic research. For 
more than three decades, IZw 18, discovered by Zwicky (1966) 
and first studied by Sargent & Searle (1970) held the record as 
the most metal-deficient galaxy known, with an oxygen abun- 
dance 12 -I- log(0/H) = 7.17 + 0.01' in its NW component and 
7.22 + 0.02 in its SE component (Izotov et al. 1999). Only re- 
cently has it been displaced by two more metal-deficient BCDs, 



' This translates into a metallicity of a; 1 /3 1 Zq with the solar oxygen 
abundance taken from Asplund et al. (2009). 



SBS 0335-052W with 12-Hlog(0/H) = 7.12±0.03 (Izotov etal. 
2005) and DDO 68 with 12 h- log (O/H) = 7.14 + 0.03 (Izotov & 
Thuan 2007). 

It was postulated that the difficulty of finding more metal- 
poor objects could be explained by thorough local mixing of 
the metals produced and released over a single starburst episode 
(Kunth & Sargent 1986). This hypothesis, implying the exis- 
tence of quasi-pristine gas beforehand, has been since chal- 
lenged by observations of uniform abundances in disconnected 
star-forming regions within a single dwarf galaxy (e.g., Skillman 
& Kennicutt 1993; Kobulnicky & Skillman 1997; Noeske et al. 
2000). Although evidence of local metal enrichment has already 
been observed (e.g., Walsh & Roy 1993), it is usually associated 
with the presence of Wolf-Rayet (WR) stars (Walsh & Roy 1989; 
Thuan et al. 1996; Lopez Sanchez et al. 2007). This hypothesis 
was later confirmed by Brinchmann et al. (2008), with the study 
of a large sample of WR galaxies from the Sloan Digital Sky 
Survey. 

Chemical abundances and metallicity are key parameters 
for the evolution of galaxies. Abundances in BCDs are usually 
probed in the ionized gas of the H n regions, as derived from op- 
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tical emission-lines. While H ii regions correspond to the present 
star-formation episode, they can also be related to previous 
episodes via the triggered star-formation mechanism, in which 
star formation take place following compression of molecular 
clouds by massive stellar winds and supernovae (e.g., Vanhala 
& Cameron 1998; Bhattal et al. 1998; Fukuda & Hanawa 2000; 
Pustilnik et al. 2001; Lee & Chen 2007; Elmegreen 201 1). H n 
regions represent however a small fraction of the galaxy in mass 
as compared to the surrounding H i region (Thuan & Martin 
1981). 

The H I region provides another way of measuring the metal- 
licity, which accounts for the bulk of the mass of a galaxy. The 
H I region sets a reference abundance pattern, caused by metal 
enrichment over long time scales and large spatial scales, that 
needs to be understood and constrained. With this in mind, the 
comparison of abundances between the H i region and the H n 
regions allows us to examine the dispersal and mixing of heavy 
elements from recent starburst episodes, and to constrain stel- 
lar evolution models. Furthermore, since future star formation 
thrives on the H i gas reservoir, the H i region makes it possible 
to study the interplay between star-formation history and metal- 
licity. While past and present star-formation episodes enrich the 
ISM in heavy elements, the latter are, in turn, able to facilitate 
and sustain star-formation through several gas cooling mecha- 
nisms (e.g., [C n] 157/im, CO, and dust radiative cooUng). The 
main questions that drive the studies of the H i region in BCDs 
are the following (1) What mechanisms control the H i region 
metallicity? (2) What is the influence of metallicity on the star- 
formation history? 

BCDs are ideal targets because they display large amounts 
of H I gas (Thuan & Martin 1981) and because the massive 
stars provide a strong far-ultraviolet (FUV) continuum on which 
absorption-lines from neutral species located along the line of 
sight are superimposed. Kunth et al. (1994) used the GHRS in- 
strument on board the Hubble Space Telescope to observe inter- 
stellar absorption lines toward the many massive stars in I Zw 18. 
The authors detected the O i line at 1302.2 A and, together with 
the neutral hydrogen content measured from the 21 cm hne, they 
estimated the metallicity of the neutral envelope to be at least 10 
times lower than that of the ionized gas in the H ii regions. Later, 
Thuan et al. (1997) obtained a GHRS spectrum for SBS0335- 
052E and found the metallicity of the neutral gas to be similar to 
that of the ionized gas. These results remained however incon- 
clusive because of possible saturation effects of the O i absorp- 
tion line (Pettini & Lipman 1995). 

A further step forward was then achieved with the launch of 
the Far Ultraviolet Spectroscopic Explorer (FUSE; Moos et al. 
2000), which allowed the observation of H i absorption-lines to- 
gether with many metallic species such as N i, O i. Si ii, P ii, Ar i, 
and Fe ii. The spectral resolution of FUSE (/l/A/l ~ 20000) and 
its sensitivity made it possible to detect numerous absorption- 
lines, some of them apparently not saturated. The neutral gas 
chemical composition was derived in several BCDs, spanning a 
wide range in ionized gas metallicity, from 1 /30 to 1 732©, with 
IZw 18 (Aloisi et al. 2003, hereafter A03; Lecavelier des Etangs 
et al. 2004, hereafter L04), SBS 0335-052 (Thuan et al. 2005), 
IZw 36 (Lebouteiller et al. 2004), Pox 36 (Lebouteiller et al. 
2009), Mark 59 (Thuan et al. 2002), NGC625 (Cannon et al. 
2005), and NGC 1705 (Heckman et al. 2001; Lee & Skillman 
2004). These investigations showed that the neutral gas of BCDs 
is not pristine; it has already been enriched with metals, with a 
minimal metallicity of ~ 1/50 Zq (see the FUSE sample analysis 
in Lebouteiller et al. 2009). The second most important result is 
that the metallicity of the neutral gas is generally lower than that 



of the ionized gas - except for the two lowest-metallicity BCDs, 
IZw 18 and SBS 0335-052E - implying that the H i region has 
been probably less processed than the H n region. 

However the interpretation of these results has been limited 
by the lack of constraints on the properties of the gas detected 
in absorption in the FUV (physical location along the line of 
sight, density structure, temperature). FUV observations select 
out dust-free hues of sight toward the massive stars, and the 
intervening gas might have peculiar properties as compared to 
the global ISM. Uncertainties also remain in deriving colunm 
densities from absorption-lines in complex systems. Most FUSE 
studies of BCDs assumed a single absorption component aris- 
ing from a single homogeneous hne of sight with the exception 
of Mark 59 (Thuan et al. 2002; see also Lebouteiller et al. 2006 
for the giant H ii region NGC 604 in M33). A more realistic 
approach would be to consider multiple lines of sight toward 
massive stars, each line of sight intersecting clouds with pos- 
sibly different physical conditions (such as turbulent and radial 
velocity) and chemical properties. Whenever individual absorp- 
tion components are unresolved, it is almost impossible to tell 
whether they are saturated, in which case the column density in- 
ferred from the global absorption line can be severely underesti- 
mated (see e.g., Lebouteiller et al. 2006). Following the sugges- 
tion of Kunth & Sargent (1986), Bowen et al. (2005) got around 
the problem by deriving the abundances in the neutral gas of the 
dwarf spiral galaxy SBS 1543+593 along a quasar line of sight. 
The authors found that the abundances agree with those in the 
ionized gas of the brightest H n region in the arms, contrasting 
with the FUSE results in BCDs. 

With the advent of the Nubble's Cosmic Origins 
Spectrograph (COS; Green et al. 201 1) we have the opportunity 
to look anew at the problem of the chemical discontinuity 
between the ionized and neutral gas in BCDs, with in particular 
the observation of the S ii multiplet, which Pettini & Lipman 
(1995) have proposed as a tracer of the neutral gas enrichment. 
The cornerstone galaxy IZw 18 is therefore once again put to the 
test, not only because of its extremely low metallicity content. 
As it turns out this object has led to discordant analysis of the 
same dataset obtained with FUSE by two independent groups of 
scientists. A03 used an individual hne fitting method and found 
abundances globally lower in the neutral gas as compared to the 
ionized gas, while L04 find no significant difference between the 
two phases. L04 find O i/H i= -4.7;^° ^, which is consistent with 
the O/H ratio observed in the H ii regions (-4.83 ± 0.03; Izotov 
et al. 1999), while A03 report a significantly different value 
with O i/H 1= -5.4 ± 0.3, a discrepancy with the former result 
that awaits clarification. This is indeed worrisome since this 
object has also a considerable weight in the study of depletion 
and that of a possible threshold for chemical enrichment. Most 
importantly, the discrepancy between the observational results 
complicates the interpretation of chemical evolution models. 
Recchi et al. (2004) found that abundances in the H i medium 
(defined by a temperature lower than 7000 K) should be similar 
to those in the H ii gas. A definitive observational constraint is 
needed to challenge the models. 

While providing accurate measurements of chemical abun- 
dances in the neutral gas, the COS observations also shed a 
new light on the physical conditions and thermal processes in 
the H I region. The absorption hne C n* /11335.7, arising from 
the fine-structure level of C^, is intimately related to the far-IR 
(FIR) coohng line [C ii] 157//m that was recently detected in 
IZw 18 with the Herschel space telescope (Lebouteiller et al. 
in prep.). We investigate in the present study the heating and 
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Fig. 1. Hubble/ ACS image of I Zw 1 8 (main body), with RBG colors de- 
fined as /, V, and B filters. Images were downloaded from the Hubble 
Legacy Archive (http://lila.stsci.edu/) and later combined. The green 
contours show the H i column density distribution from Lelli et al. 
(2012) with 2" resolution (beam size in the bottom left). Contours are 
drawn for 3 (dashed), 6, 9, 12, and 15 x 10"' cm"^. The COS observa- 
tion, represented by the black circle (2.5" diameter), is centered on the 
NW star cluster. 



Table 1. COS datasets. 



# 


a 


6 


Grating/f" 




(J2000) 


(J2000) 




11523 


09:34:02.00 


+55:14:28.0 


G130M/3.4, G160M/7.3, 








G185M/8.1 


12028 


09:34:02.00 


+55:14:28.0 


G130M/8.2, G160M/11.4 


11579 


09:34:01.97 


+55:14:28.1 


G130M/15.5 



Notes. The exact coordinates for observations 1 1523 and 12028 are un- 
known because a peak-up was performed beforehand to target the UV 
brightest region. Coordinates are J2000. 

Exposure time in units of xlO'' s. 
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Fig. 2. Processing of the observed spectrum (a) consists of correcting 
for reddening and for Lya absorption in the Milky Way and in the high- 
velocity cloud (b); fitting and removing hya from IZw 18 (c); and nor- 
malizing the spectrum (d). The residual emission at Lyo- wavelength is 
due to geocoronal H i. 



cooling mechanisms in the H i region in order to understand the 
present and future star-formation history of the galaxy. 

Our study is organized as follows. We present the observa- 
tions in Sect. 2. Column densities are determined in Sect. 3. After 
the description of the nebular models in Sect. 4, we examine the 
physical conditions in the H i region in Sect. 5. We then derive 
and interpret elementary abundances in Sect. 6. Finally, we dis- 
cuss chemical evolution scenarios in Sect. 7. 

2. Observations and data processing 

2.1. Datasets 

To first approximation, IZw 18 is composed of two young, mas- 
sive star clusters, NW and SE, which ionize the surrounding 
medium. We used COS on the Hubble Space Telescope (HST) 
to observe the brightest region in the FUV, NW, in April 2010 
(program 1 1523) and later in December 2010 (program 12028). 
Both observations were made with the circular aperture, with a 
diameter D ~ 2.5", corresponding to ~ 230 pc at the distance 



of IZw 18 (~ 19Mpc; Aloisi et al. 2007). Observations were 
centered on NW (see Fig. 1 and the observation log in Table 1). 
Target acquisition consisted of a slew to the specified coordi- 
nates of NW followed by a peak-up in ultraviolet light. Although 
both observing programs made use of the FUV gratings G130M 
and G160M, the two sets of exposures cover somewhat differ- 
ent wavelength range. Together, they span ~ 1135 - 1760 A. 
We also use the archival dataset from program 1 1579 which ob- 
served with the G130M grating (~ 1 132 - 1433A). No peak-up 
was performed for the latter observation and the telescope was 
slewed to the requested coordinates of NW. The G185M spec- 
trum was not used because of the low data quality^. 



- The G185M spectrum contains in particular the Si n /11808.0 line, 
with an oscillator strength about 50 times lower than the weakest Si ii 
line in our study, /il546.7. However, the low S/N of the spectrum, 
due to the almost fully illuminated aperture (Sect. 3.2) and to cross- 
contamination between the spectral stripes of the grating on the detec- 
tor, does not allow setting a useful constraint on the Si ii column density. 
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The pipeline reduction procedure is continually being im- 
proved, and we obtained our spectra well after the observations 
were made (July 2011) from the Hubble Data Archive, which 
makes an "on-the-fly" reduction using the best calibrations avail- 
able. We combined exposures from a given program to form a 
single, co-added spectrum. We did not merge the spectra from 
the three programs because the pointing was shghtly different. 
This results in differences in the spectra, with for instance the 
C IV /11548, /11550 doublet, which is formed in stellar winds and 
in the H ii region. 

Among the 3 available datasets, the 11579 observation has 
the longest exposure time below ~ 1433A, while the 12028 ob- 
servation has the longest exposure time above (Table 1). The 
12028 observations should be superior to the 11523 observa- 
tions by virtue of their longer exposure times. Indeed, the higher 
signal-to-noise ratio (S/N) in the 12028 spectra enabled us to dis- 
cern broad, shallow stellar absorption lines that we had not de- 
tected before in the other datasets. Also, the 12028 observations 
were accomplished in a single 5-orbit visit, so we are assured 
of a constant telescope pointing throughout the visit, whereas in 
program 11523, the G130M exposures and G160M exposures 
were obtained in different visits. 

In some exposures, geocoronal emission (in particular from 
the 1 1302, 1305, 1306A multiplet) is contaminating the inter- 
stellar line profiles. We ignored the corresponding exposures for 
the wavelength range around the geocoronal emission. 

As a preparation for the fitting of the interstellar lines, we 
prepared each COS spectrum in several steps as illustrated in 
Fig. 2. In the four steps, we: 

- Corrected the observed spectrum for foreground Galactic ex- 
tinction by dust, using the color excess E{B - V) - 0.032 
(Schlegel et al. 1998). This correction is used for the nor- 
malization by the stellar model (see below); 

- Corrected for Lya absorption in the galaxy and for Lyor ab- 
sorption by the high-velocity cloud in the line of sight to 
IZw 18, following Kunth et al. (1994); 

- Measured the column density of Lya in I Zw 1 8 and removed 
Lyor (Sect. 3.4); 

- Normalized by the stellar spectrum (continuum and absorp- 
tion in photospheres; Sect. 2.2). Figure 3 shows the spectra 
along with the normalized stellar model. Final spectra are 
shown in the Appendix. 

2.2. Contamination by stellar absorption 

The FUV spectrum of IZw 18 is clearly composite, with con- 
tributions from stars, H n region, and H i region. Absorption in 
stellar photospheres needs to be removed since they may con- 
taminate interstellar absorption-lines. We calculated model spec- 
tra of young stellar populations using grids of high-resolution 
(0.2 A) non-LTE photospheric spectra (900-2000 A) by Lanz & 
Hubeny (2003, 2007) and Geneva evolutionary tracks (Lejeune 
&Schaerer2001). 

High-resolution imaging studies (Hunter & Thronson 1995; 
Heap et al. 2013) indicate that IZw 18 is undergoing continu- 
ous star formation rather than an instantaneous burst. Also, it 
is clear from deep HST/ACS observations that IZw 18 has a 
widely extended, older stellar population (Aloisi et al. 2007; 
Contreras Ramos et al. 2011). However, it is equally clear that 
young, massive stars are the dominant contributors to the UV 
flux of the NW component. Observations of individual stars 
in IZw 18 with WFC2 (Hunter & Thronson 1995) and STIS 
(Brown et al. 2002; Heap et al. 2013) on Hubble indicate that 



Table 2. Absorption systems. 



System 


species 


v" 

(kms-') 


(kms-') 


IZw 18 


neutraF 


~ 755 


~ 17.8 


Milky Way 


neutraP 


« -10 


« 22 




neutraF 


« -90 


«7 


Milky Way 


ionized (S m. Si m, P m) 


* -60 


* 85 


A VI cloud'' 


neuttal*^ 


* -171 


*21 



^"^ Radial velocity. Turbulent velocity. The species in the neutral 
component are those Usted in Table 5. High- velocity cloud. 



the NW massive star cluster has a range of ages. The bright- 
est stars in the UV are also the youngest, and their inferred 
masses range up to 150 M©. We therefore assumed continuous 
star formation and a Salpeter initial mass function with lower 
and upper mass limits, 1 - 150 Mq, to calculate the spectrum 
of the stellar population for ages (i.e., the duration of star for- 
mation up until now) ranging from to lOOMyr. We then con- 
volved each calculated spectrum with a rotational broadening 
function (vsin; = 150kms"') and with the instrumental func- 
tion (Sect. 3.2). From a comparison of these models to pho- 
tospheric lines in the COS spectrum, we derive a continuous- 
star-formation age of 10 - 15Myr and a stellar metalUcity 
log(Z/Zo) ~ -1.7. 

Using the age and metallicity derived above, we calculated 
the overall spectral energy distribution (SED) of stars in the NW 
component using the grids of low-resolution, LTE spectra by 
Castelli & Kurucz (2003). This grid and the TLUSTY/synspec 
grid only cover effective temperatures up to 50000 K, somewhat 
lower than the highest temperatures of massive stars on the zero- 
age main sequence (65000 K), and definitely lower than highly 
evolved stars such as WR stars. It is therefore likely that our 
model SED underestimates the ionizing flux of the stellar popu- 
lation. 

Figure 3 shows the stellar models overplotted on each 
dataset. We find that the contamination from stellar fines is neg- 
ligible for most of the interstellar lines from IZw 18 and also 
from the local absorption systems (Sect. 3.1). One exception is 
the contamination of S ii /11250.6 in the Milky Way by the stellar 
C m line. The stellar model was used to normalize the observed 
spectrum so as to isolate the purely interstellar absorption. 

3. Column density determination 

3.1. Method 

Once the spectrum has been corrected for stellar absorption 
(Sect. 2.2), several interstellar absorption line systems can be 
distinguished along the line of sight to the massive stars in 
I Zw 18 (Table 2). The main interstellar lines are listed in Table 3. 
The atomic data, in particular the oscillator strength values, 
were taken from Morton (2003), except for P ii and S ii for 
which we used Federman et al. (2007) and Podobedova et al. 
(2009) respectively. The I Zw 1 8 component is visible at a ve- 
locity of » 755 km s"^ Two local components are required at 
« -lOkms"^ and « -90km s"^ In addition, a component at 
~ -60kms"' is detected that corresponds to ionized species in 
the Milky Way. Finally, the high-velocity cloud "A vi" in com- 
plex Ais detected at ^ -171 kms"' (Hulsbosch & Wakker 1988; 
Wakker & van Woerden 1997). 

We noticed a relative shift in the velocity of the neutral com- 
ponent within IZw 18 between the datasets: 756 ± 3kms"' for 
12028, 750 ± 2kms-i for 11523, and 768 ± 3kms-i for 11579 
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Fig. 3. COS spectra of IZw 18. The Lya absorption (from IZw 18, the Milky Way, and the high- velocity cloud) was removed and the flux was 
normalized to the local continuum. The numerous absorption features often come in pairs: the weaker - redshifted - absorption feature is from 
IZw 18; the stronger feature arises in the Galaxy. A non-LTE, photospheric model spectrum (no winds) of the stellar population (lOMyr age with 
a continuous star formation hypothesis) and log(Z/Zo) = -1.7 is shown in blue. While many features such as C iii AllTS and the Si iv /11400 
doublet appear to be well matched by the stellar model, other high-ionization features such as He ii /1I640 and the C iv /il549 doublet are not well 
matched. The origin of these lines (e.g., WR stars, H ii region...) is unclear. 



(Table 4). The same relative shifts are observed for the Milky 
Way component, suggesting an instrumental origin. We attribute 
this difference to a different distribution of the UV flux in and 
outside the aperture. As a comparison to our velocity determi- 
nations, the Ha observation of the NW component, probing the 
ionized gas in the H ii region, is 742 + 7kms"'(Petrosian et al. 



1997) while the H i 21 cm line has a velocity of x 740 km s"' 
(Lelli et al. 2012). Table 4 also shows the radial velocities in- 
ferred from each species individually. 

The column densities were calculated with the line fitting 
routine Owens (Lemoine et al. 2002). Owens was originally de- 
veloped to solve the complexity of many interleaved absorption 
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systems in the FUSE spectra. Although line blending is not a 
significant issue in the COS spectram of I Zw 1 8, we used Owens 
because of its ability to perform a simultaneous fit of the absorp- 
tion line profiles from any number of absorption systems. Each 
gaseous system is defined by (1) a radial velocity Vrad, (2) a tem- 
perature T, (3) a turbulent velocity b, and (4) a set of column 
densities from selected species. In practice, the species within 
a given absorption system are assumed to share the same set of 
parameters (vrad, T,b). 

We assume in the following that the absorption arises from 
a single broad component as opposed to a collection of narrow 
components. Within this hypothesis, the column density deter- 
mined from saturated lines highly depends on the b parameter. 
Systematic uncertainties on b include errors on the instrumen- 
tal resolution (Sect. 3.2) or the fact that species in the neutral gas 
might not share exactly the same b value. If we choose to neglect 
these systematic uncertainties, the simultaneous fit of lines with 
a wide range of optical depths provides constraints on both b and 
the column density. We verified that the additional presence of 
saturated lines in the simultaneous fit does not bias the column 
density determination if the column density is otherwise already 
well constrained by weak lines. 

Other effects might also be at work that introduce system- 
atic errors on the column density. For example, the velocity 
distribution could include a number of individual narrow ab- 
sorption components that we caimot resolve. The integrated, 
unresolved, line profile may appear optically thin even though 
narrow components are saturated (the so-called "hidden satura- 
tion"). Although the absorption spectra of metal-poor objects, 
with their low metal column densities, is in principle less af- 
fected than more metal-rich sources, the hidden saturation effect 
is non-linear and can result in underestimating the total column 
density by several factors (e.g., Lebouteiller 2005; Lebouteiller 
et al. 2006). Observationally, this effect can be minimized by in- 
cluding lines with a low oscillator strength / (e.g., C ii*, S n, 
P II in our spectra). Another method is to verify that, for a given 
species, lines with a wide range of / values can be fitted by a 
single column density and that the apparent optical depth scales 
with / (e.g, N I, S ii, Fe n in our spectra). Nevertheless, the col- 
umn density of some other species was determined only from 
strong lines. We note in particular C ii /11334.5 (optical depth 
in the Une center tq =s 76), O i .11302.2 (tq ~ 34), and all the 
Si II lines, with the weakest being /11526.7 (tq ~ 9). For these 
species, we rely on the hypothesis that the absorption does not 
arise from a collection of narrow components, and we keep in 
mind that the quoted errors could be underestimated. 

The absorption lines from the various systems are blended 
with each other throughout the COS range. The following lines 
in IZw 18 are contaminated by local absorption systems, Fe ii 
Al 142.4, Si II il 190.4, /II 193.3, Mn ii il 199.4, and S ii .11250.6. 
The Owens algorithm allows us to use these lines to determine 
column densities under certain conditions, either other lines 
from the same species are isolated and/or the column density 
of the local component is well constrained by other lines. In 
the FUV spectrum of IZw 18, Fe ii. Si ii, Mn ii, and S ii have 
other isolated lines and the local absorption systems are well 
constrained by a variety of lines throughout the COS range. 

In order to derive column densities, Owens performs a 
minimization over a number of spectral windows containing one 
or several Unes. While Owens enables the parametrization of 
the local continuum around the lines, the normalization we per- 
formed (Sect. 2.1) results in a flat continuum that does not need 
to be further parametrized. For each specie, we included every 
line in the COS wavelength domain. This includes the lines that 



Table 3. Main interstellar Unes used in this study. 



Line 


(A) 


/ 




Line 


A 

(A) 




/ 


c 


I* 


1657.907 


4 


71e-02 


S I 


1425 


030 


1.25e-01 


C 


I* 


1657.379 


3 


56e-02 


S I* 


1323 


516 


3.06e-02 


c 




1656.928 


1 


49e-01 


S I 


1316 


543 


3.45e-02 


c 


I* 


1656.267 


5 


88e-02 


S I* 


1310 


194 


2.89e-02 


c 


I* 


1560.682 


9 


59e-03 


S I* 


1277 


216 


4.93e-02 


c 




1560.309 


1 


28e-02 


S I 


1270 


780 


5.51e-02 


c 


I* 


1329.123 


1 


60e-02 


S I* 


1253 


325 


2.90e-02 


c 


I* 


1329.100 


2 


60e-02 


S I 


1247 


160 


3.24e-02 


c 


I* 


1329.085 


2 


13e-02 


S II 


1259 


518 


1.66e-02 


c 




1328.833 


7 


58e-02 


S II 


1253 


805 


1.21e-02 


c 




1280.135 


2 


63e-02 


S II 


1250 


578 


5.43e-03 


c 


I* 


1279.891 


1 


43e-02 


CI I 


1347 


240 


1.53e-01 


c 


I* 


1277.513 


2 


lOe-02 


CI I 


1335 


726 


3.13e-02 


c 


I* 


1277.283 


7 


05e-02 


CI I 


1188 


774 


7.01e-02 


c 




1277.245 


8 


53e-02 


Mn n 


1201 


118 


1.21e-01 


c 


I* 


1261.122 


1 


47e-02 


Mnn 


1199 


391 


1.69e-01 


c 


I* 


1260.926 


1 


35e-02 


Mn n 


1197 


184 


2.17e-01 


c 




1260.735 


5 


07e-02 


Mn II 


1162 


015 


1.02e-02 


c 




1193.030 


4 


09e-02 


Fe n 


1611 


200 


1.38e-03 


c 


I* 


1193.008 


2 


91e-02 


Fe n* 


1618 


468 


2.20e-02 


c 




1157.910 


2 


12e-02 


Fe n 


1608 


451 


5.77e-02 


c 


I* 


1157.770 


1 


83e-02 


Fe n* 


1267 


422 


1.69e-02 


c 


n* 


1335.708 


1 


15e-01 


Fe n* 


1266 


677 


8.99e-03 


c 


n* 


1335.663 


1 


28e-02 


Fen 


1260 


533 


2.40e-02 


c 


n 


1334.532 


1 


28e-01 


Fe n* 


1148 


277 


8.28e-02 


N 




1200.710 


4 


32e-02 


Fe n 


1144 


938 


8.30e-02 


N 




1200.223 


8 


69e-02 


Fe n 


1143 


226 


1.92e-02 


N 




1199.550 


1 


32e-01 


Fe n 


1142 


366 


4.01e-03 


N 




1134.980 


4 


16e-02 


Fe n 


1133 


665 


4.72e-03 


N 




1134.415 


2 


87e-02 


Nin 


1773 


949 


6.22e-03 


N 




1134.165 


1 


46e-02 


Niii* 


1754 


813 


1.59e-02 





I* 


1304.858 


4 


78e-02 


Nin 


1751 


916 


2.77e-02 


o 




1355.598 


1 


16e-06 


Nin 


1709 


600 


6.88e-02 


o 




1302.168 


4 


80e-02 


Ni n* 


1788 


491 


2.52e-02 


Al 


n 


1670.787 


1 


74e+00 


Ni n* 


1748 


289 


1.40e-01 


Si 




1526.707 


1 


33e-01 


Nin 


1741 


553 


4.27e-02 


Si 




1260.422 


1 


18e+00 


Nin 


1703 


412 


6.00e-03 


Si 




1193.290 


5 


82e-01 


Nin 


1502 


148 


1.33e-02 


Si 




1190.416 


2 


92e-01 


Nin 


1467 


756 


9.90e-03 


P 




1787.648 


6 


04e-02 


Niii* 


1500 


434 


4.77e-02 


P 




1782.829 


1 


13e-01 


Nin 


1467 


259 


6.30e-03 


P 




1774.949 


1 


69e-01 


Nin 


1454 


842 


3.23e-02 


P 




1679.697 


5 


36e-02 


Nin 


1415 


720 


3.19e-03 


P 




1674.595 


3 


25e-02 


Nin 


1412 


866 


3.55e-03 


P 




1381.476 


3 


16e-01 


Niii* 


1411 


065 


6.28e-02 


P 




1379.428 


2 


19e-01 


Nin 


1381 


685 


5.08e-03 


P 




1377.073 


1 


lOe-01 


Niii* 


1381 


286 


1.15e-01 


P 




1152.818 


2 


72e-01 


Nin 


1370 


132 


7.69e-02 


S 




1473.994 


7 


30e-02 


Nin 


1345 


878 


7.69e-03 


s 


I* 


1433.309 


5 


85e-02 


Niii* 


1335 


201 


1.63e-01 


s 


I* 


1433.278 


1 


69e-01 


Nin 


1317 


217 


1.46e-01 


s 




1425.188 


3 


65e-02 


Nin 


1308 


866 


5.41e-03 



Notes. X is the rest wavelength, / is the transition oscillator strength. 
The * denotes the transitions from the fine-structure level. 



are not detected, since the latter help providing a reUable upper 
limit to the column density. The fits were performed on the 3 
observation sets separately because of possible systematic dif- 
ferences in the spectra. 
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Table 4. Radial velocities (2cr error bars). 



Species 


11523 


12028 


11579 


Global 


750 ±2 


756 ±3 


768 ±3 



Ci 

Ci* 

Cn 

Cn* 

Ni 

Oi 

Aln 

Sin 

Sin* 

Pi 

Pn 

Si 

Sn 

Cli 

Mn n 

Fe n 

Nin 



755+5 


763+« 


769+5 


139*1 


749+^' 


748+J 




737:^ 


769!' 


761!!^ 






748::^ 


753^6 




15TI 


763:^5 


765+^ 



747:11 



749+10 



72r;« 

759:^ 

743::f 



763:« 



765!« 
770!^ 

769:;^ 



Notes. The first row gives the radial velocity of the neutral component, 
as derived from the simultaneous fit all the species in the table (as- 
suming that the species share the same radial and turbulent velocities; 
Sect. 3.1). For each species, we also provide the radial velocity calcu- 
lated by isolating the species in its own component, i.e., the radial ve- 
locity and the turbulent velocity are not constrained by other species. 



3.2. Line broadening 

The observed line width is a combination of the intrinsic broad- 
ening (thermal and turbulent) with the instrumental line spread 
function. The latter is a critical parameter upon which the b value 
determination relies. Systematic errors on the line spread func- 
tion could lead to wrong diagnostics on the line saturation and, 
therefore, to unreliable column density determinations. 

In the case of a point-source well-centered in the aper- 
ture, the instrumental full width at half maximum (FWHM) is 
~ 0.08A. If the source has extended emission or if it is a col- 
lection of point-sources, the spectral resolution is significantly 
degraded^. Depending on the source distribution in the aperture, 
the Une profiles might be shifted and/or skewed, but we found 
the lines to be Gaussian and symmetric, suggesting a rather uni- 
formly distributed light in the aperture. 

We estimated the effective instrumental spectral resolution in 
two ways. First, we measured the one-dimensional flux profile 
of stars in the NW component as derived from a FUV image ob- 
tained by Brown et al. (2002). The profile looks roughly gaussian 
with a FWHM= 1.25", i.e., half the size of the COS aperture. 
The flat-topped profile of geocoronal Lya emission, which uni- 
formly fills the COS aperture, has a FWHM= 0.84A, so the ex- 
pected instrumental profile is expected to have a FWHM= 0.4A. 
The second estimate of the spectral resolution makes use of 
the data itself. For this, we considered the sharpest Unes in the 
COS spectrum and also the line blends for which the 2 lines are 
barely separated (in particular at x 1134A and a; 1260A). These 
profiles give a reliable constraint on the instrumental FWHM, 
« 0.45 - 0.50A. We use 0.50A (« QOkms"') for tiie fine fitting 
of the 11523 and 12028 datasets (Sect. 3.1). A shghtly larger 



COS instrument handbook, http://www.stsci.edu/hst/cos. Sect. 5.9. 



FWHM was required for the 1 1579 dataset, 0.55 A, which could 
be due to pointing jitter over the long exposure time. For com- 
parison, the point-Uke source FWHM with the FUSE telescope 
is 20kms-i (Hebrard et al. 2002; Wood et al. 2002). 

We used the instrumental resolution above to calculate the 
intrinsic line broadening. For the latter, we could not distinguish 
between the thermal and turbulent components. We assumed 
that turbulence is the dominant source of broadening by setting 
the temperature to 1000 K, corresponding to a thermal broaden- 
ing of < Ikms"'. The turbulent velocity value we determine 
for IZw 18 is » 17.8 + 3kms~' for the 12028 dataset, and 
18.2 ± 8kms~' for the 11579 dataset, in good agreement with 
the value 17.3 ± 5kms"' found by L04. The 11523 dataset did 
not converge on a well constrained b value, which is likely due 
to the relatively low S/N ratio of this observation and also to the 
presence of spurious features. We forced the b value in the 1 1523 
dataset to be « 17.8 ± 3 km s"' . 

3.3. Heavy elements 

The fitting of the absorption line profiles is shown in the 
Appendix. Final column densities are given in Table 5. The error 
bars were calculated by using the Ax^ method, which gives the 
variation oix^ ~Af^> where is the lowest value corre- 
sponding to the global minimum over the entire parameter space. 
Since we use the information from all available lines in the wave- 
length range, the error bar can be significantly smaller than when 
considering a single fine. We refer to L04 and Lebouteiller et al. 
(2006, 2009) for more details. The comparison between the col- 
umn densities we derive and the FUSE results asks for caution. 
The large FUSE aperture (30" x 30") includes the entire main 
body of IZw 18 while COS only observed the NW component. 
However, most of the UV-bright stars are located toward NW 
(Brown et al. 2002), and the corresponding Unes of sight pre- 
sumably dominated the FUSE spectral continuum. Therefore, 
we do not expect large differences between the column densities 
derived with COS and with FUSE. In the following we provide 
individual coimnents for each species. 

The C n column density is not well constrained because the 
only C II line, /11334.4, is heavily saturated. Our determination 
relies on the confidence that the b value is well constrained by 
other lines from other species included in the fitting and on the 
fact that hidden saturation effects are negligible (Sect. 3.1). The 
column density of C n* is much better determined. C i and C i* 
lines could be present in our spectra but they are either barely 
detected or not detected. The contribution from C i and C i* to 
the total C colimm density is negligible. 

The N I /II 134 triplet, observed in the 11579 dataset, is not 
saturated (see also L04), while the /11200 triplet is close to satu- 
ration. The N I column density determination from the 3 datasets 
are in good agreement. They fall somewhat above the value from 
L04 and agree well with that from A03. 

O I is observed in the /II 302.2 line (oscillator strength / = 
0.048) in the 11523 and 11579 datasets and in the /11355.6 fine 
(oscillator strength / = 10"^) in all datasets. The former line 
is contaminated by geocoronal O i emission (Sect. 2.1), and, al- 
though the contaminated exposures were ignored, we cannot ex- 
clude that a weak contamination remains in the other exposures. 
The /11355.6 line is never detected and provides a large upper 
limit to the O i column density. L04 and A03 measured different 
O I column densities from the same spectrum, mostly because of 
the set of O i lines used for the fitting. L04 ignored the /11039 
line because of contamination by a terrestrial airglow and be- 
cause of the line saturation. Our O i column density determina- 
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Table 5. Column densities. 



11523 



12028 



COS 
11579 



Adopted 



L04 



FUSE 

A03 (±lo-) 



Model" 



Hi 



21.34 



21.34 ±0.04 21.34 



21.34 ±0.04 



21.30 ±0.10 21.35 ±0.10 



21.34 



12-85-- 

< 13.18 
(15.38tU:^«) 

1 ^ CQ+0.14 
^■^•O-^_0.07 

(14.45ti^) 

< 18.08 

(13.25!«:^?) 
(14-53t«) 

< 11.69 

< 12.51 

< 12.86 

< 12.95 

1 a 9Q+0.14 
^■^•^■^-0.27 

(14.49!«;^«) 

1 1. 4O+0.08 



Ci 

Ci* 

Cii 

Cii* 

Ni 

Oi 

Oi* 

Aln 

Sin 

Sin* 

Pi 

Pii 

Si 

Sn 

Cli 

Mn n 

Fe n 

Nin 



< 12.91 
(15-46t«) 

(14 55+'' '>2-) 
(16.05!«:^^) 
(< 13.31) 
(13.32:»i^) 
(14.76t») 

< 12.31 

< 12.10 
12 91+0 21 
12 93+o '8 
14 75+0 09 



13 10+0 05 
(14 71+0 081 
n 37+0- 11 

1 J.J ,Q 



< 13.42 

< 12.97 
(15.37:0-3^) 
13.76!«;» 
(14.59!0:0^) 

(16.00!°:^) 

(< 13.87) 



(14-5ltg) 

< 12.06 

< 12.18 

12 49+0 31 

< 12.30 
(14.56t08)<i 

< 13.05 

< 13.30 

(i4.v6!!J:«?) 

1 ■% 47+0.04 



< 13.40 

< 13.00 
15.35 ± 0.30 
13.80 ±0.10 
14.55 ±0.10 
* 16.37* 

(< 13.90) 
13.32 ±0.30 
14.52 ± 0.30 

< 11.90 

< 12.50 
12.92 ± 0.30 

< 12.30 
14.73 ±0.10 

< 13.05 
13.10 ±0.50 
14.75 ± 0.20 
13.42 ±0.10 



14 79+0.14 

16 60"^°-^° 



14.80!«:| 



14.60t|^ 



14.44 ± 0.04 
15.98 ± 0.26 



14.81 ±0.07 



< -13.60 



15.09 ± 0.06 



15.49 (-0.05) 

13.53 (-0.37) 

14.54 (-0.01) 



13.32 (-0.19) 
14.86 (-0.00) 



14.74 (-0.00) 



13 


12 (- 


-0 


00) 


14 


66 (- 


-0 


00) 


13 


33 (- 


-0 


01) 



Notes. Error bar.s are Icr error bars unless noted otherwise. Note that the error bars for the 1 1523 dataset are likely underestimated because of the 
tight constraint on b for this dataset (Sect. 3.2). The values in parentheses could be affected by systematic uncertainties on the h value or due to 
spurious features (see text). The O i column density in dataset 12028 is an upper limit since it is determined using the undetected /11355 line. 
Column density calculated by the Cloudspec model (Sect. 4). The value in parentheses gives the correction when the emission component is 
accounted for (Sect. 4.2). Estimated using S n column density as a ttacer (Sect. 6.2). '•'^^ No coverage. Noisy profile. 



tion is closer to that of A03, but we consider our value as highly 
uncertain because of saturation. 

Aluminum has a single strong line in the COS range at 
1670. 8A. The line is unfortunately saturated, and the Al n col- 
umn density determination relies on the constraints on the b 
value. We adopt an average column density from the datasets 
11523 and 12028. 

All the Si n lines are saturated. The Si n column determina- 
tion is thus relatively uncertain since it depends strongly on the b 
value determination and on the assumption that absorption arises 
from a single broad component (Sect. 3.1). We use both 12028 
and 11579 datasets for the Si ii column density determination, 
and ignore the 1 1523 dataset because of its relatively lower S/N. 
Our adopted column density is somewhat lower than the FUSE 
measurements by L04 and A03. We thus keep in mind that our 
Si n column density might be underestimated and adopt a large 
error bar that reconciles our value with the FUSE values. 

The strongest P n line (Al 152.8) is barely detected. Our col- 
umn density determination falls well below the upper limit from 
A03. The strongest P i line at 1381.5A is not detected, and the 
contribution from P i to the total P column density is negUgible. 

S II is observed in the A1250.8, /11253.8, and /11259.5 lines. 
The main constraint on the S ii column density is given by the 
weakest line, /11253.8. It is well detected in the 12028 dataset, 
while the profile is more noisy in the 2 other datasets, espe- 
cially in 11579. We adopt an average column density between 
the 11523 and 12028 datasets. 

Chlorine is expected as CI i and CI ii in the neutral phase. 
CI II has no lines in the COS range and it was not detected in the 
FUSE spectra. CI i has numerous lines below ~ 1400A with the 
strongest at 1347.24A. The latter is however not detected in any 
of the datasets. 



Manganese is observed via several Mn ii lines with strong os- 
cillator strengths. Unfortunately all the lines are heavily blended 
with saturated (or close to being saturated) N i hues from I Zw 18 
and from the Milky Way. Therefore, we consider that the Mn ii 
column density is highly uncertain and may be affected by sys- 
tematic uncertainties. 

Iron is observed via several Fe n lines spanning a wide range 
of oscillator strengths. The presence of weak Unes should re- 
sult in a reliable Fe n column density determination. Our col- 
umn density determination falls between the value from L04 and 
A03. It agrees slightly better with L04 while it is not compatible 
within errors with the value from A03. The difference between 
L04 and A03 is due to the lower b value (< 8 km s~') assumed 
by A03 for the Fe n lines. 

Nickel is observed via many weak Ni ii lines, most of which 
are not detected. We stress that the column density determina- 
tion of Ni n is better than 5cr because it uses the information 
from all the lines we included in the fitting. The column density 
determination from the 3 datasets are in remarkable agreement. 

3.4. Hydrogen 

The column density of H i is an important quantity since metal 
abundances are tied to it. However, it is a difficult parameter to 
measure correctly. The blue wing of the Lyor profile is contami- 
nated by numerous metal lines and by Lya absorption from the 
Galaxy and from a high-velocity cloud along the line of sight 
to IZw 18 (Sect. 3.1). Furthermore, the red wing of Lya is con- 
taminated by nebular continuum emission (2-photon and Balmer 
continuum) arising in the H ii region of IZw 18 and by the N v 
/11238, /11240 wind lines from the ionizing stars in IZw 18. 
Finally, the continuum level is uncertain, because Lya is such 



8 



Lebouteiller et al.: Chemical enrichment and physical conditions in I Zw 18 



Table 6. H i column density determinations. 



Value 


Method 


Instrument" 


Ref. 


Absorption 


21.34 ± 0.04 


Lya (NW) 


HST/COS (~ 2.5") 


(1) 


« 21.32 


LyyS 


FUSE (30" X 30") 


(2) 


21.54 ±0.06 


Lya (NW) 


lUE & HST/GHRS (2") 


(3) 


21.34 ±0.07 


Lyman series 


FUSE (30" X 30") 


(4) 


21.32 ±0.09 


Lyman series 


FUSE (30" x30") 


(5) 


22.3 


Lyii' (peak) 


HST/STIS (52" X 0.5") 


(6) 


Emission 


« 21.2 


21 cm (NW) 


VLA (8") 


(7) 


^ 21.20 


21 em 


VLA (20") 


(8) 


21.53 ±0.15 


21 cm (NW) 


VLA (5") 


(8) 


Si 21.48 


21 cm (NW) 


VLA (2") 


(8) 



References. (1) This study; (2) Vidal-Madjar et al. (2000); (3) Kunth 
et al. (1994); (4) A03; (5) L04; (6) Brown et al. (2002); (7) van Zee 
et al. (1998); (8) LeUi et al. (2012). 

The value in parentheses indicate either the aperture (absorption 
measurements) or the beam size (emission measurements). 



a broad line. In principle, both red and blue wings can be con- 
taminated by stellar Lya absorption, but this contamination turns 
out to be insignificant (see also Lebouteiller et al. 2008). 

We worked with the COS spectrum processed through step 
3 (i.e., corrected for Lyes' in the Milky Way and in the high- 
velocity cloud; Sect. 2.1). We measured the Lya column den- 
sity in IZwl8 by multiplying the observed spectrum by e^, 
with T = A^CH i)xlO->'-3'' X (A - /to)-2 (Jenkins 1971). The 
column density that best canceled the observed absorption is 
A^(H i)=2.2x 10^' cm"^. We find a similar value for the 3 datasets. 
We list in Table 6 the various H i column density determinations 
in the Uterature, either derived in absorption from the Lyman se- 
ries or in emission using the hyperfine- structure 21 cm line. Note 
that the comparison between observations in absorption and in 
emission should account for the simple geometrical effect that 
the absorption only probes the gas in the foreground of the UV- 
bright stars. 

The H I spatial distribution observed by Lelli et al. (2012) is 
particularly informative. The high-spatial resolution (2" beam) 
reveals that the NW region lies in a H i hole (see also Fig. 1). 
The NW region is located between the southeast H i peak, with 
A^(H i)« 1.5 X 10^^ cm~-^, and a region to the north with a rather 
flat column density distribution (» 6 X 10^' cm~^). Interpolating 
between these 2 values, we would expect an H i column density 
of ~ 10^^ cm"^ toward NW, which is significantly larger than our 
observed value. As discussed in Lelli et al. (2012), the young 
stellar cluster in NW likely created an H i hole, either due to 
ionization or to a blow out of the ISM. Assuming the H i hole is 
due to ionization, the column density of ionized hydrogen should 
be on first approximation 10^^ - 2.2 X 10^' = 7.8 x 10^' cm"^. 
We discuss further the distribution of ionized gas in Sect. 4. 

4. Nebular models 

4.1. Geometry and physical conditions 

Models were constructed to compute the ionization structure of 
the nebula and the relative emission and absorption strength of 
the intervening medium along the line of sight to UV bright stars. 
We made the foUowing simplifying assumptions: (1) the ioniz- 
ing stars of NW can be treated as a single, central source; (2) 
nebular emission arises over the whole volume of H ii region, but 
nebular absorption arises only in material along the line of sight 



to the central source; and (3) there is a clear-cut boundary be- 
tween the H II and H i regions. The latter assumption is problem- 
atic. The photoionization model of IZw 18 by Pequignot (2008) 
indicates that the covering fraction of the H ii region is only 40%, 
so the bulk of the photoionizing radiation must escape into the 
H I region, and deep observations of IZw 18 show Ha emission 
in arcs outside the star-forming region and in a diffuse halo 3 kpc 
across (Dufour & Hester 1990). 

Spectral analysis of the H i cloud of IZw 18 is much more 
comphcated than for normal (i.e., star-less) damped Lya systems 
(DLAs). For example, the spectral profile of the C n 1334-5 res- 
onance doublet is formed by absorption by relatively lower-mass 
and/or cooler young stars, and by emission from the outer parts 
of the H II region, as well as by absorption in the H i cloud. 
Even Lya, nominally formed in the H i cloud is not immune to 
contamination. Balmer continuum emission and 2-photon emis- 
sion from the H ii region affect the long-wavelength side of the 
line, and absorption by H" atoms in the atmospheres of cool 
stars can affect both the short- and long-wavelength wings of 
Lya (Sect. 3.4). 

We used the software package, cloudspec (Hubeny et al. 
2000), which combines the photoionization code Cloudy 
(Ferland et al. 1998) with the spectral-synthesis code synspec 
(Hubeny & Lanz 1995) to compute the emission and absorp- 
tion spectrum of an object. The physical model is a medium 
photoionized by a central source and the cosmic background. 
Cloudy computes the ionization structure and emission spec- 
trum (both Unes and continuum) of the medium. The synspec 
part of cloudspec then calculates the level populations and the 
attenuation spectrum produced by such a medium along the fine 
of sight towards the central ionizing source. 

We ran cloudspec through the H ii region and into the H i 
region until the H i column density reached its observed value 
toward the NW component, 2.2 x 10^^ cm"^ (Sect. 3.4). Figure4 
shows the radial profile of the electron temperature and density 
and the sharp change in the physical conditions between the H n 
and the H i region. Close in to the ionizing stars, the nebula is 
quite hot, due to the low abundance of oxygen, the main el- 
ement for coUisional cooling in the ionized gas. The electron 
temperature decreases outward from the ionizing stars with an 
abrupt drop near the edge of the H ii region, and outside the 
Stromgren sphere, the temperature reaches down to lOOK. The 
electron density drops by 3 orders of magnitude at the edge of 
the Stromgren sphere, because most of the free electrons sup- 
plied by hydrogen are no longer available. 

The NW component is dominated by the H ii region which 
sits an H i hole either due to hydrogen ionization or to a blow 
out of the ISM (LeUi et al. 2012). In our model, we find that 
the column density of ionized hydrogen should be logA^(H i)= 
21 .79, i.e., almost a factor 3 larger than the observed H i column 
density. Following the hypothesis that the total (ionized-nneutral) 
hydrogen column density is flat across the main body of IZw 18, 
we would expect the ionized hydrogen column density toward 
NW to be around log A^(H i)~ 21.89 in order to compensate for 
the relatively lower H i column density (Sect. 3.4). Our models 
are in good agreement with this value and support the idea that 
the gas was not blown out but mainly ionized. 

In assuming spherical symmetry, the models are highly sim- 
plistic. The ionizing stars in IZw 18 are offset from Ha emis- 
sion shell (Hunter & Thronson 1995; Brown et al. 2002), and 
the H n region has a knot of high density. Nevertheless, the mod- 
els help us to understand important physical processes operating 
in IZw 18 (Sect. 5) and to guide the abundance analysis of the 
H I region (Sect. 6). 
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Fig. 5. Cloudspec model spectra. Top - Stellar and nebular (H i and H 11 region) component normalized to the stellar continuum (blue) and 
emission component (red). Both component spectra have been convolved with a gaussian line-spread function having a FWHM= 0.5A to simulate 
the observations. Major absorption lines formed in the H i region are labeled on top. Most other strong lines, e.g. the Si iv /1 1400 doublet, C iv 
/il549 doublet, are higher ionization lines formed in the H 11 region. The absorption component of Lya completely erases the Lyo- emission formed 
in the H 11 region. Bottom - Simple model spectrum for the case CF x AF = 0.1 (top) and for CF x AF = 0.4 (bottom), with CF the covering 
fraction and AF the aperture fraction (see text). The dotted line is the absorption spectrum, and the bold line is the net spectrum. The component 
and net spectra have been convolved with a gaussian line-spread function having a FWHM= 0.5A. 



4.2. Spectral line formation in the H 11 region 

Singly ionized ions are formed by recombination of higher 
ionization stages inside the H n region near the edge of the 
Stromgren sphere. Recombination and subsequent cascade down 
to the ground state in the H 11 region creates emission that works 
to reduce the apparent strength of absorption lines arising in the 
H I region. 

Figure 5a shows the emission and absorption components of 
the model line spectrum as derived by cloudspec. The observed 
spectrum is the net. As can be seen from the figure, the emis- 



sion is strong relative to the absorption for C 11* /11335.7. How 
important is contamination by emission lines? The answer de- 
pends on two factors: the covering fraction of the H 11 region, 
CF, and the fraction of the H 11 emission region subtended by 
the entrance aperture to the spectrograph, AF. Pequignot (2008) 
estimates the covering fraction of the H 11 region in IZw 18 at 
CF - 0.4, based on the strength of Hfi emission to Lyman con- 
tinuum flux of the model best fitting the observations. The COS 
aperture fraction is not known, but is certainly much less than 
1 . Figure 5b compares a segment of the model UV spectrum for 
two cases, CFxAF - 0.1 (top) and CFxAF - 0.4 (bottom). It is 
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Fig. 4. Electron temperature (top) and density (bottom) profile of the 
simple model approximation of IZw 18. The abscissa is the radial dis- 
tance from the ionizing stars in cm. 



clear that the latter case disagrees with observation by generally 
overestimating the strength of emission, making the ground-state 
absorption lines too weak, and the fine-structure lines appearing 
in emission (particularly C ii* and Si ii*). The agreement is much 
better for the case where CF x AF = 0. 1 . In our abundance anal- 
ysis, we have therefore neglected the effects of emission except 
for the C n doublet at 1334-5 A. 

Contamination by emission cannot always be neglected. In 
fact, the spectra of star-forming galaxies at high-redshift often 
show emission in the fine- structure lines (Pettini et al. 2000; 
Shapley et al. 2003). This emission is easily identifiable because 
these galaxies have outflows. The H i region is expanding with 
respect to the central, ionizing stars, so the absorption compo- 
nent, which arises along the line of sight to the stars, is blue- 
shifted. The emission component is formed over the outermost 
shell within the H ii region, so it is centered on the stellar ra- 
dial velocity. Expansion of the H ii region simply broadens the 
line. Thus, the red half of the emission reaches the observer un- 
scathed, but the blue half of the emission is at least partially nul- 
lified by the blue-shifted absorption component. The net profile 
looks like a P Cygni profile. The emission spectrum of I Zw 1 8 is 
as real as in high-redshift, star-forming galaxies. The difference 
is that there is no detectable velocity gradient in IZw 18, so the 
emission is not separable from the absorption component. 

While the emission component in the H n region is negligi- 
ble for most lines, absorption can also arise from the H ii region. 
A few resonant lines originating in the ionized gas are seen in 
absorption in the FUV spectra of IZw 18 and provide an inde- 
pendent, direct, estimate of the ionization structure of the nebula. 
For instance, A03 measured the Fe m /II 122 line in the FUSE 
spectrum of IZw 18 and set an upper limit to the fraction of Fe ii 
in the ionized gas of 5%. A03 also cautioned that a large fraction 
of Si n absorption might arise in the ionized gas. We are now able 
to test this hypothesis, with the Si in /11206.5 line, which we de- 
tect at a velocity of ~ 755 + 50 km s ' in the COS spectra. This 
line partly originates from stellar atmospheres and partly from 
the H II region. The observed profile of the Si m line is com- 
plex and it depends on the stellar contribution that was removed 
(Sect. 2.2). The column density determination also depends on 
the b value for which we have little information in the ionized 
gas. Assuming that the line is not saturated, we find a consis- 
tent Si m column density in the 3 datasets with ~ 2 x 10'^ cm"-^. 
What is the fraction of the Si n column density we measured 
(3.3 X 10'^ cm"^) arising in the ionized gas? In the warm ionized 
medium in the Milky Way, about 90% of sihcon is expected into 



Si II and 10% into Si iii (Sembach et al. 2000). Assuming this 
result is apphcable to IZw 18, we find that only ~ 40% of the 
measured Si n column density might arise from the neutral gas. 
However, the situation is likely different in IZw 18 where sili- 
con in the ionized gas is expected mostly into higher ionization 
stages (Sect. 6. 1 ; see also A03). We assume in the following that 
most of Si n should lie in the H i region gas. 

Unfortunately, no S ra line could be used to perform a simi- 
lar analysis with sulfur The brightest S iii line, Al 190, is blended 
with Si II. Nevertheless, we consider that the ionization correc- 
tion smaller than for Si and that most of S n Ues in the H i region. 



5. Physical conditions and thermal equilibrium in 
the H I region 

The C n* column density represents only » 3% of the total car- 
bon column density (Table 5), yet C ii* is key to understanding 
the physical conditions, in particular the thermal processes to- 
ward the NW region. The ground state of C"^ is split into the 
fine-structure levels ^Pi/2 and ^P3/2. The absorption line C n 
/11334.5 originates from the ^Pi/2 level while C n* /11335.7 orig- 
inates from ^P3/2. The upper level is populated by collisions with 
H*^, H2, and e~. The spontaneous emission between the two fine- 
structure levels is observed in the far-infrared emission line [C 11] 
157//m and dominates over collisional de-excitation as long as 
the density is lower than the transition critical density. 



5.1. Electron fraction in ttie H i region 

The population of the C^ fine-structure level constrains the elec- 
tron density toward the NW region. Following Jenkins et al. 
(2000a), the volume-density ratio «(C^*)/n(C^) provides the 
electron density as long as collisions with electrons dominate the 
population of the fine-structure level. Collisions with H2 should 
be unimportant in IZw 18 given the lack of diffuse H2 (Vidal- 
Madjar et al. 2000). Furthermore, the collision strength with e~ 
is larger than the collision strength with H atoms by a factor 
of ^ 150 for T = 6000 K. Colhsions with will therefore 
dominate as long as ne/«H remains larger than ~ 0.75% (e.g., 
Lehner et al. 2004). Potential sources of free electrons are soft 
X-rays, cosmic rays, and ionization of C". We approximate the 
number density ratio n(C^*)/n(C^) with the column density ra- 
tio Rqii = N(C^*)/N(C'^) and, by assuming that electrons are 
the dominant collision partners, we derive an upper limit on the 
electron density 



121 



cn 



721 (e) 2e-^«=/^-/?, 



cm 



(1) 



cii 



where A21 is the radiative decay probability from the fine- 
structure level from the ^P3/2 level to the ^Pi/2, T'exc is the ^P3/2 
level excitation temperature (« 91 K), T is the electron temper- 
ature and 721(e) is the collisional deexcitation rate coefficient, 
which we take from Goldsmith et al. (2012). 

For a temperature range of [1000 - 8000] K, we find < 
[10 - 20] /?cncm"^. Using the observed coluiim density ratio 
Ren = 2.8 X 10"^, we calculate He < 0.4 cm"^. Since the C 11 
absorption-line is saturated, we also checked that the result holds 
when using S n as a replacement for C 11. The S n lines are not 
prone to saturation, and both S n and C n are the dominant ion- 
ization stages in the H i region, so we expect the ionization frac- 
tions of S into S"^ and C into C"^ to be well correlated over a 
wide range of conditions. Carbon is more readily depleted on 
dust grains than sulfur (Sofia et al. 1997, 1998), but rather than 
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assuming a depletion factor, we assume that the abundance ratio 
in the neutral gas equates that in the ionized gas (C/S)Hn ~ 0.98 
(Pequignot 2008). The latter assumption also implies that C and 
S disperse and mix in the ISM over similar timescales (Sect. 6.4). 
We can then approximate Rqii with 



1 



CII 



jV(C+*) 

N(S^) (C/SW 



(2) 



from which we derive a somewhat lower upper limit of tie ^ 
0.2cm-3. 

On the one hand, the upper limit on the electron density 
we determine is significantly lower than the value measured in 
the ionized gas of the star-forming regions derived in Dufour & 
Hester (1990) and in the Ha arc derived by Izotov et al. (2001a), 
on the order of lOcm"^. This suggests that C ii* arises in the 
diffuse ionized gas or in the H i region. On the other hand, the 
upper limit on He is significantly larger than the electron density 
in the H i region inferred from our Cloudy models (~ 0.01 cm"^; 
Sect. 4). The models suggest that the C ii* emission component 
can be neglected and thus does not impact our estimate of Rca 
(Sect. 4.2). However, it is plausible that collisions partners other 
than e~ contribute to populating the C"^ fine-structure level, in 
particular H atoms. We can use the same method as Eq. 1 to de- 
rive an upper Umit on the hydrogen density, i.e., assuming that 
hydrogen atoms are the dominant collision partners. For a tem- 
perature range of [100 - 1000] K, we find hh < [30 - 40] cm ^ 

If we assume that the H i region has a density around or 
lower than 10 cm"^, we must conclude that the electron fraction 
in the H i region is Ukely significant, possibly reaching a few per- 
cents. Such values would lie many orders of magnitude above the 
fraction of electrons provided by ionization of C" alone (equal 
to the carbon abundance relative to hydrogen), suggesting that 
the H I region is partially ionized, with electrons created by UV 
photons, soft X-rays, or cosmic rays originating from the star- 
forming regions. We investigate further the electron fraction in 
Sect. 5.2. 

5.2. Theoretical cooling rate from the H i envelope 

In order to examine further the physical conditions in the H i 
region, we now compare the observed cooling rate to the theo- 
retical rate expected from the H i envelope, with H atoms and 
free electrons as collisions partners. The observed cooling rate 
per H atom is calculated following Pottasch et al. (1979) and 
Gryetal. (1992), with 

L = A2rh y2i ^ - 2.89 x 10"^° erg s'^ (H'^) (3) 

where V2i is the transition frequency and where we have again 
approximated the number density ratio with the column density 
ratio. We find Lcos ~ 8 x 10"^^ erg s"' (H"^). We can rewrite the 
cooling rate per surface unit of the galaxy. 



L' ^ 2.89 X 10"^" NiC^*) erg s"' cm 



-2 



(4) 



which leads to L'^q^ 



1.8 X lO^^'erg s"' cm"'. Assuming the 
source is extended, the intensity per surface unit of the galaxy 
is the same as the intensity per surface unit on the detector, so 
we can directly relate L' to the corresponding 157 jum emission 
intensity /(157/im) in erg s"' cm"^ sr"^ via: 



1(157 fim) = L' 



4n' 



where <1> is the fraction of the beam filled with [C ii] emis- 
sion. Assuming O = 1, we find /cos(157jum) ~ 1.5 x 



10 ^ erg s ^ cm ^ sr" 



1 = 1.5 X 10-1" Wm-2 sr-i. This is the ob- 



served cooling rate toward the NW component. 

For the theoretical cooling rate in the H i envelope, we 
consider a medium with most free electrons produced by cos- 
mic rays, soft X-rays, or escaping UV photons (see Sect. 5.1). 
Following Crawford et al. (1985) and Madden et al. (1993), the 
cooling rate expressed in W m"-^ sr'^ can be written as 



/th(157//m) = 5.8 X 10 



,-7 



(Zq) (1021 cm-2) 



Fc+ O y. 



(6) 



where Z is the abundance of carbon relative to the solar value 
(C/H= 2.69 X 10-''; Asplund et al. 2009), A^h is the hydrogen 
nuclei column density, Fq* is the fraction of carbon into ions 
along the line of sight, <t> is the fraction of the beam filled by C"^ 
gas, and y is the dimensionless factor defined as 



■1 ' 



(7) 



(5) 



2e-9i/J- 

y ~ 

1-1-2 -I- ine/na,e + nH/ncr,H)' 

where Wcr.e and Wcr.H are the critical densities for collisions with 
electrons and hydrogen (which we take from Goldsmith et al. 
2012). We assume in the following Fq* = 1 and we define x as 
the electron fraction x = rie/nii. Figure 6 provides /th(157/im) as 
a function of density for various temperatures and electron frac- 
tions. When the ionization fraction is large, proton impact could 
contribute to the ^Ps/i level excitation, but here we neglect this 
contribution since we focus on the H i region. It can be seen that 
/th(157yum) flattens for densities > lOcm"-' for large x values 
(medium significantly ionized), and > 10^ cm"^ for low x values 
(purely neutral medium), which is a consequence of the criti- 
cal densities for colUsions with electrons and hydrogen atoms. 
/th(157//m) depends little on temperature for T greater than a 
few hundreds K. 

In IZw 18, the abundance of carbon relative to solar is Z = 
10-2, where we have used sulfur abundance as a substitute and 
assumed the C/S ratio in the ionized gas (see Sect. 5.1). We first 
consider a mostly neutral medium, with T ~ 1000 K and use 
the H I column density derived in Sect. 3.4 (IQ-^^-^'^ cm"2) From 
Fig. 6, we calculate that the observed cooling rate /cos(157/im) 
can be reproduced by a wide range of conditions, from purely 
neutral with nn ~ 15 cm-^, to a medium with a significant elec- 
tron fraction with njj « 5cm-^ and x « 2%. Assuming a ho- 
mogenous medium with a low clumping factor, such number 
densities imply small H i scales on the order of > 100 pc, which 
is compatible with an H i shell, left over after the NW star cluster 
blew out and ionized the pre-existing cloud (Sect. 4.2). 

For larger electron fractions, the medium is diffuse and ion- 
ized and we need to account for the column density of H^. 
Following the hypothesis that the total hydrogen colutmi density 
toward NW is « IO22 cm-2 (gg^.^ 3 4-)^ gjj^ ^^^^ ^jjg cooUng 
rate can be reproduced by an electron density of k 0.04 cm~^. 
However, such a low number density would result in an unreal- 
istic nebula size (> 50kpc). 

In summary, although the lack of knowledge on the hydrogen 
number density in the H i region prevents us from deriving a 
precise electron fraction, our results are compatible with most 
of the C II* absorption arising in the H i region with likely a 
significant electron fraction. 

How does the cooling rate derived from C n* compare 
with the predicted value from the Cloudy models? The mod- 
els (Sect. 4) predict that radiative de-excitation of the fine- 
structure level dominates the cooling in the H i region, with 
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x= -4 -3 -2 -1 




"h [cm-'] 



where Esfr is the SFR per unit area and [Z] represents the 
D/G normalized to the MW value. While [Z] should in prin- 
ciple scale with the gas metallicity, we use instead the D/G as 
measured with the Herschel Space Telescope (2000;^ [3*] times 
lower than the MW value; Remy et al. in prep). We assume 
SsfrCMW) ~ 4x 10-3 Mq yi--' kpc^^ (Rana 1991) and L(MW) * 
2 X lO^^^'ergs"' (H)"' (Lehner et al. 2004), and we finally ob- 
tain Esfr(IZw18) ~ 0.3 ± O.2M0yr"' kpc"2. As an alterna- 
tive measurement of the SFR, we use the empirical relation 
between the gas surface density Egas and Ssfr (e.g., Kennicutt 
1998). The H i column density is used as a tracer of Sg^s, 
with a coiTection factor of ~ 4 due to the ionized gas along 
the line of sight (Sect. 4.1). We find Igas ~ 87.6 Mopc"^ and 
Ssfr(IZw18) ~ 0.13+;}20MQyr-ikpc°2. Both these determi- 
nations are close to the values found by Aloisi et al. (1999), 
with 10^2 - 10"' Moyr"' kpc"^ and Petrosian et al. (1997), 
with 0. 14 Mq yr"' kpc"^ (coiTected to a distance of 19 Mpc). The 
slightly higher SFR found from the C 11* cooling rate, if genuine, 
might be due to an underestimated D/G because of the presence 
of metal-free gas on the line of sight (Sect. 7.1) and/or to a larger 
photoelectric efficiency in IZw 18 as compared to the MW (see 
also Sect. 5.4). 



Fig. 6. The 157 yum cooling line intensity is calculated from Eq. 6 and is 
plotted as a function of the hydrogen number density and of the elec- 
tron fraction (x = rig/nu)- Note that the result should be scaled by the 
fraction of carbon in C^. The result should also be scaled by the carbon 
abundance with respect to the solar value and by the column density of 
hydrogen nuclei in units of 10^' cm"^. For temperatures 50 K and 100 K, 
only the curve for x = is shown. 

■^'modei ~ 3.5 X 10"^^ erg s~' (H"'), which is a factor of ^ 2 lower 
than the observed value Lcos ~ 8 x 10"^^ erg s"' (H"'). While 
we assumed a constant hydrogen density of 10 cm"^ in the mod- 
els, a lower density would result in a higher temperature, with 
other coolants contributing significantly to the total cooling rate, 
such as [O i] 63 fim or [Si ii] 35 fj.m. This would suggest that the 
total cooling rate in the models is likely underestimated by sev- 
eral factors. Since our models assume a thermal equilibrium, the 
heating rate is also underestimated by the same amount. We dis- 
cuss in Sect. 5.4 the heating mechanisms and how to reconcile 
the models with the observed C ii* cooling rate. 

5.3. Cooling rate and SFR 

The observed cooling rate derived from Eq. 3 is Lcos ~ 8 x 
10"^^ erg s ' (H '). This value is significantly lower than in the 
diffiise medium of our galaxy (a: 10"^^ - 10"^^ erg s"' (H"'); 
Lehner et al. 2004) or in other nearby galaxies (e.g.. Madden 
et al. 1993), but it is close to the average cooling rate in DLAs 
(a 10^2^ erg S-' (H-i); Wolfe et al. 2003a). Assuming that the 
photoelectric effect on dust is the dominant gas heating mech- 
anism (see Sect. 5.4), the cooling rate scales with the dust-to- 
gas ratio (D/G), the photoelectric efficiency, and the FUV field. 
Since the latter is a good gauge of the star-formation rate (SFR) 
per unit physical area, the similar cooling rate between IZw 18 
and DLAs could in principle imply a similar surface SFR den- 
sity (Wolfe et al. 2003b). Following Bowen et al. (2005), we can 
estimate the SFR from C ii* by using the Milky Way (MW) as a 
reference 

2sfr(IZw 18) ^ [Z]-' Ssfr(MW), (8) 



5.4. Heating meclianisms 

In the warm ionized medium, photoionization of H should dom- 
inate the heating (see Pequignot 2008 for IZw 18), with likely 
extra heating from photoelectric effect on dust and from dissipa- 
tion of interstellar turbulence (Reynolds et al. 1999). The heat- 
ing mechanisms in the H i region are comparatively much less 
constrained. The gas heating rate in the ISM of our Galaxy is 
dominated by photoelectrons ejected from polycyclic aromatic 
hydrocarbons (PAHs) and small dust grains (Wolfire et al. 1995; 
Weingartner & Draine 2001), but low metallicity environments 
could show fundamental differences, with for instance the low 
PAH abundance of PAHs (e.g., Wu et al. 2006). 

Ultraviolet photons escaping the H ii region play an impor- 
tant role in the H i gas heating through photoionization and 
through the photoelectric effect. Dufour & Hester (1990) discov- 
ered a faint halo of Ha emission extending at least 2.7 kpc from 
the star-forming massive star clusters of IZw 18. Furthermore, 
Pequignot (2008) found that only 30 - 40% of ionizing pho- 
tons from the NW and SE massive star clusters is absorbed by 
surrounding H ii regions. The leakage of ionizing photons into 
the neutral medium should result in a heating effect, affecting the 
relative population of the ground and fine-structure levels of C ii, 
as well as the strength of the C ii* /11335 fine-structure line. 

The total heating rate predicted by the Cloudy models in the 
H I region (Sect. 4) is dominated by cosmic rays, and amounts 
to r ~ 3.5 X 10"^^ erg s"' cm"3. Photoelectric effect heating 
is comparatively negligible, with ~ 2% of the total heating 
rate'*. However, our simple model has important limitations. 
First, the heating rate is likely underestimated by several factors 
(Sect. 5.3), which could be due an underestimated photoelectric 
effect heating rate. The latter is highly dependent on the model 
geometry and parameters since it is proportional to the incident 
UV radiation field, the photoelectric efficiency, and the D/G ra- 
tio (e.g.. Bakes & Tielens 1994; Weingartner & Draine 2001). 
Apart from the remarkably low D/G in IZw 18 (Remy et al. in 



^ If the D/G is simply scaled to the metallicity of IZw 18, the photo- 
electric effect heating rate becomes the dominant heating mechanism, 
contributing to half of the total heating rate, itself increased by a factor 
of 2. 
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prep.), the low photoelectric heating rate predicted by our mod- 
els is due to the relatively weak radiation field in the H i region 
(Go < 100, expressed in Habing units; Habing 1968), and to the 
low photoelectric efficiency, due to the large grain charging pa- 



rameter (Go 



Iq5 ]^i/2 cjjj3-) jjjg photoelectric effect 



might become the dominant heating mechanism in the models if 
a more realistic approach was taken. We list below some of the 
most important parameters, 

- The radiation field in the H i region depends on the ISM 
topology and in particular on the covering factor of the H ii 
region. The UV photons likely escape the star-forming re- 
gion to reach and heat the gas on large spatial scales across 
the galaxy. 

- The electron density in the models might be underestimated 
(Sect. 5.1). A larger electron density will increase the photo- 
electric efficiency, resulting in a larger heating rate. 

- The grain size distribution could be biased toward the small- 
est sizes in low-metallicity dwarf galaxies (e.g., Galliano 
et al. 2005), with likely a significant impact in the gas heating 
rate and efficiency. 

- The D/G ratio is uncertain, between 700 and 3000 times 
lower than the MW value. 

The study of the far-lR [C ii] line and its relation to the IR lumi- 
nosity will be important to confirm the importance of dust in the 
gas heating process (Lebouteiller et al. 2013). 

Since our model appears to underestimate the heating rate 
(Sect. 5.3), we explored further the heating of the H i gas through 
the use of the "extra heating" parameter in Cloudy. The ex- 
tra heating represents all heating sources that are not well con- 
strained or not modeled (e.g., escaping ionizing photons, com- 
pressional heating, dissipation of turbulence, mechanical energy 
transfer from stellar winds and supernovas...). We assumed no 
radial dependence for this parameter. Figure? shows that the 
electron temperature, initially at lOOK in absence of extra heat- 
ing, jumps from 300 K to 7000 K when the extra heating rate is 
increased from 10"^^ to 10"^^ ergs"' cm"-'. The higher temper- 
ature brings increased collisional excitation of to the fine- 
structure level, and hence, stronger absorption of C ii* /11335.7. 
As expected, the strengths of the emission components to the C ii 
doublet are hardly affected. The emission appears very strong - 
enough to overwhelm the absorption component of C ii /11334.5, 
but the COS aperture only subtended a small fraction of the H n 
region, so the emission components are significantly weaker (see 
Sect.4.2). Nevertheless, the unknown strength of the emission 
component adds to the uncertainty on the heating rate as derived 
from the observed C ii* profile. 

Assuming an extra heating rate of 
we are able to reconcile our models with the cooling rate de- 
rived from C II* (Sect. 5.3). The electron temperature Tp in the 
H I region with such extra heating increases from 100 K to 
1000 K. The electron temperature is unfortunately difficult to 
measure observationally, and, although the column density ratio 
A^(Si^*)/A^(C^*) provides in principle the gas kinetic temperature 
(Howk et al. 2005), our upper limit on the Si ii* column density 
is too large to provide any constraint. 

By including the extra heating, the main IR cooling lines 
predicted by the models are [C ii] (~ 40% of the total cooling 
rate), [O i] 63 yum (~ 30%), and [Si ii] 35 fim (~ 26%). We 
wish to recall here that these values correspond to the diffuse 
medium, and that far-lR emission-line intensities might include 
other components (e.g., photodissociation regions). The O i* 
/11304.9 absorption-line, arising from the ^^Pi fine-structure level 
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Fig. 7. Effect of extra heating in the H i region. The electron tempera- 
ture (lop panel) and the fraction of C ii and C ii* to total carbon (middle 
panel) are shown as a function of radius to the ionizing stars. The pre- 
dicted line strength is shown in the bottom panel. The models assume 
unity for both the covering fraction of the H n region and the fraction of 
the emission subtended by the COS aperture (see Sect. 4.2). 

of O", is unfortunately undetected. Following Eq. 3, we calcu- 
late an upper limit on the [O i] cooling of < 10"^^ erg s"' (H"'), 
which is significantly larger than the cooling by [C ii] inferred 
from C II* (Sect. 5.3). Therefore, we cannot examine in detail 
the effect of extra heating on the relative cooling line contribu- 
tions. At best, the large upper limit on the cooling by [O i] is 
compatible with the extra heating. 



6. Chemical abundances 

6.1. Dominant ionization stages in the H i region 

In order to calculate the abundances relative to atomic hydrogen, 
all the ionization stages in the neutral phase should be accounted 
for. Fortunately, the COS range gives access to the lines of the 
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dominant ionization stages in the neutral gas, namely C n, N i, 
O I, Mg n, Si n, P ii, S ii, Mn ii, Fe n, and Ni n. Only chlorine 
is expected about equally between CI i and CI n. Most studies 
of DLAs assume that singly ionized metals such as C, Al, Si, 
S, Fe, and Ni reside exclusively in the H i region so that ele- 
mental abundances can be obtained from the ratio of the ion col- 
umn density to the H i column density (e.g., Pettini et al. 2000). 
However, this assumption is not always correct for H i clouds 
with embedded star-forming regions. To illustrate, we have con- 
structed a simple Cloudy model of the gaseous medium sur- 
rounding the ionizing stars (Sect. 4). The model is simplistic in 
having a constant total hydrogen density, 10 cm"^, but in other 
respects - metallicity, ionizing radiation field, star-formation 
rate, and H i column density - the model is appropriate for 
IZwl8. 

Figure 8 compares the model spectral energy distribution of 
the ionizing stars incident on the surrounding gaseous medium 
and after passage through the H ii region. No H-, N-, or O- 
ionizing radiation gets through to the H i region, so oxygen and 
nitrogen are neutral in the H i region. However, stellar radiation 
longward of the H i Lyman Umit, although diluted by distance 
from ionizing stars, has basically the same relative flux distribu- 
tion except for the addition of very strong Lya and other, much 
weaker nebular features. Hence, other elements such as C, Al, 
Ni, Mg, Fe, Si, S, and CI will be ionized in the H i region. 

From Fig. 9, we can see that in the H i region, the dominant 
ionization state of metals with ionization edges longward of the 
H I Lyman limit is singly ionized. We use the column densities 
of the corresponding species to determine the total abundance. 
A fraction of these species also arises in the H ii region where it 
can lead to emission (recombination of higher stages; Sect. 4.2) 
and absorption. 

A caveat should be that partially-ionized regions might ex- 
ist on the lines of sight that our model cannot account for. The 
existence of such regions might affect the ionization structure of 
species with a large photoionization cross-sections (PiCS). Even 
if the ionization potential is greater than 13.6 eV, the PiCS can 
be large enough for ionization to occur while hydrogen remains 
atomic. This effect is unimportant for thick clouds well shielded 
from the FUV radiation, but becomes significant for thin clouds 
or clouds embedded in a diffuse UV field. Among the species 
observed in absorption in galaxies, Ar i is the species the most 
affected because of its large PiCS (e.g., Sofia & Jenkins 1998; 
Jenkins 2013). More relevant to our study, N i and O i also have 
a PiCS larger than that of H i (Verner et al. 1996), but also lower 
than that of Ar i by a factor of > 3, so we expect a much smaller 
correction than for Ar. Most importantly, the ionization fraction 
of N and O are coupled to the ionization fraction of H via reso- 
nant charge exchange reactions (see also Jenkins et al. 2000b). 
For all the other species^, the first ionization potential is lower 
than 13.6 eV, and the second ionization potential is well above, 
so we always expect them as single-charged ions in the neutral 
phase even in partially-ionized regions. 

In the following, we compare our abundances to the solar 
values from Asplund et al. (2009) except for argon for which 
we choose Lodders (2003). For the H ii region abundances, we 
use the results of Pequignot (2008) and adopt error bars that en- 
compass the determinations by Izotov et al. (1999). Following 
LebouteiUer et al. (2004), we define the abundance difference 
between the H n and H i region as: 

Sm(X) = log(X/H)„„ - log(X/H)Hi, (9) 

' Except Mn n, which is the only species for which we could not find 
the PiCS. 



Table 8. Comparison between abundances in the H i and H n region. 



X 


(X/H)e« 




[X/H]hi 


[X/Hlnn 






SmOQ 


C 


-3.57 ± 


0.05 


-2.41 ± 0.30 


-1.88 ± 


0. 


16 


0.53 ± 0.34 


N 


-4.17 ± 


0.05 


-2.62 ±0.12 


-2.23 ± 


0. 


,07 


0.39 ±0.14 


o 


-3.31 ± 


0.03 


-1.66 ±0.15* 


-1.48 ± 


0. 


06 


0.18 ±0.16 


Ne 


-4.07 ± 


0.10 




-1.54 ± 


0. 


11 




Al 


-5.55 ± 


0.03 


-2.47 ± 0.30 










Si 


-4.49 ± 


0.03 


-2.33 ± 0.30 


-1.62 ± 


0. 


,30 


0.71 ±0.42 


P 


-6.59 ± 


0.03 


-1.83 ±0.30 










S 


-4.88 ± 


0.03 


-1.73 ±0.10 


-1.55 ± 


0. 


,06 


0.18 ±0.11 


CI 


-6.50 ± 


0.30 


< -1.79 










Ar 


-5.38 ± 


0.08 


-2.42 ± 0.17" 


-1.65 ± 


0. 


,13 


0.77 ± 0.21 


Mn 


-6.57 ± 


0.04 


-1.77 ±0.50 










Fe 


-4.50 ± 


0.04 


-2.09 ± 0.20 


-1.85 ± 


0. 


,20 


0.24 ± 0.28 


Ni 


-5.78 ± 


0.04 


-2.14 ±0.11 











Abundances in the present-day solar photosphere from Asplund 
et al. (2009) except for Ar for which we take Lodders (2003). 
Calculated from S n column density and assuming S/0=-1.64 
(ionized gas value from Pequignot 2008). Argon abundance is taken 
from A03 and L04, with 12 + log(Ar/H) = 4.20 ± 0.15 (Table 5). 



where (X/H)hii and (X/H)hi are the abundances of the element 
X in the H n region and H i region respectively. We also define 
the abundance relative to the solar value as: 

[X/Y] = log(X/Y) - log(X/Y)o, (10) 

where (X/Y)o is the solar abundance reference. 

Figure 10 shows the abundances in the neutral gas (see 
Table 7 and Sect. 6) as compared those in the ionized gas. The 
line fits are shown in the appendix, including the line profiles as 
they would appear if the H i region abundances would be similar 
to those in the H ii region. 

L04 and A03 examined the same FUSE spectrum and the dif- 
ferences are due to a different method of measuring the column 
densities and also to the choice of lines (Sect. 3.3). Our method 
is similar to L04's method, i.e., a simultaneous fit of the lines. 
The differences in the abundances between our study and L04 
and A03 directly reflect the differences in the metal colunrn den- 
sity determinations since otir H i column density is similar to 
their's and since no ionization correction was used in our study 
nor in L04 and A03. 

6.2. 0/H and the a-elements 

a-elements are produced in short-lived high-mass stars and are 
often used to trace the metallicity of galaxies. The majority of 
oxygen is produced by He and Ne buming during the hydro- 
static burning phase with negligible contribution from the explo- 
sive stages. Silicon is produced by O burning (hydrostatic and 
explosive) and Ne burning, while sulfur is produced by O burn- 
ing. Since the a-elements are made in the same massive stars 
that make oxygen, the abundance ratio between two a-elements 
should depend little on metallicity (see Izotov et al. 1999). 

The sulfur abundance is tightly correlated with that of oxy- 
gen in BCDs (Izotov et al. 1999). The S/O ratio in the ionized gas 
of IZw 18 is identical flie solar ratio, witii [S/0]= -0.07 ± 0.08. 
Within the error bars, we find that S is barely underabundant in 
the neutral phase with (5hi(S) 0.18 + 0.11. Figure A.2 shows 
that the line profiles are not consistent with the H ii region abun- 
dance. Sulfur is not expected to be depleted on dust (except 
in molecular clouds) and ionization corrections are negligible 
(Sect. 6.1), so we consider that this result implies a genuine un- 
derabundance of sulfur in the neutral gas. In other words, the 
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Fig. 8. Cloudy model spectral energy distribution of 1 Zw 18. The flux distribution of the stellar population incident on the inner edge of the H n 
region is shown in white, the stellar flux transmitted to the H i region in gray. The wavelength of the ionization edge of relevant species is given 
along the top of the figure. The H ii region absorbs virtually all the H i Lyman continuum and converts it to H i Lycy, nebular continuum, and line 
emission (all shown in red). 



Table 7. Chemical abundances. 



X 


This study 


L04 (FUSE) 


A03 (FUSE) 


H II re 


»ion 


Adopted (X/H)hii 


c 


6.02 + 0.30 






6.30 + 


0.10 (3), 6.55 (5), 6.54 + 0.10 (6), 6.47 + 0.15 (7) 


6.55 + 0.15 


N 


5.21+0.11 


4.92 + 0.20 


5.09 + 0.11 


5.58 + 


0.60 (1), 5.61 + 0.03 (4), 5.60 (5) 


5.60 + 0.05 





7.05 ±0.15" 


'■-'"-0.60 


6.63 + 0.28 


7.13 + 


0.05 (1), 7.17 + 0.03 (2), 7.22 + 0.01 (4), 7.21 (5) 


7.21+0.05 


Ne 






6.38 + 


0.02 (1,2), 6.45 + 0.02 (4), 6.39 (5) 


6.39 + 0.05 


Al 


3.98 + 0.30 












Si 


5.18 + 0.30 


5.50+0.30 


5.46 + 0.12 


5.61 + 


0.23 (3), 5.89 (5), 5.65 + 0.22 (7) 


5.89 + 0.30 


P 


3.58 + 0.30 




< 4.25 








S 


5.39 + 0.11 






5.57 + 


0.40 (1), 5.75 + 0.06 (2), 5.59 + 0.05 (4), 5.57 (5) 


5.57 + 0.05 


CI 


< 3.71 












Ar 




4.15 + 0.30 


4.25 + 0.13 


4.87 + 


0.15 (1), 5.10 + 0.04 (2), 4.92 + 0.02 (4), 4.97 (5) 


4.97 + 0.10 


Mn 


3.76 + 0.50 












Fe 


5.41+0.20 


5.30+0.15 


5.74 + 0.12 


5.95 + 


0.90 (1), 5.41 + 0.09 (4), 5.63 (5) 


5.63 + 0.20 


Ni 


4.08 + 0.11 













References. (1) (09 34 02.1 +55 14 25 J2000) Thuan & Izotov (2005); (2) (NW) Izotov et al. (1999); (3) (NW) Izotov & Thuan (1999); (4) (09 34 
02.2 +55 14 25.06) Izotov et al. (1997), (5) Pequignot (2008), (6) Gamett et al. (1997), (7) Garnett et al. (1995). 

Notes. Abundances are given as 12 + log (X/H). 

Calculated from S ii column density and assuming S/0=-1.64 (ionized gas value from Pequignot 2008). 



metallicity of the neutral gas is between 51% and 85% that of 
the H II region. 

The silicon abundance determination is fairly uncertain in 
the H II region, with the values from Pequignot (2008) and 
Izotov et al. (1999) differing by almost a factor 2. Taking at face 
the value from Pequignot (2008), Si/O in the ionized gas is in 
good agreement with the solar ratio, with [Si/0]= -0.14 ± 0.33. 
Considering the similar nucleosynthesis process of these ele- 



ments, this indicates little or no depletion of Si on dust grains 
in the H ii region. It must be noted that Peimbert & Peimbert 
(2010) used the relatively smaller silicon abundance found by 
Garnett et al. (1995) to conclude that ~ 39% of Si might be 
locked in dust in IZw 18. The silicon abundance in the H i re- 
gion is uncertain because of the Si ii line saturation; our deter- 
mination is somewhat lower than the column density derived by 
both L04 and A03 (Table 5). Within errors, we find that silicon is 
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Depth [cm] 



Depth [cm] 



Fig. 9. Ionization structure of the simple model approximation of IZw 18 NW. The abscissa is the radial distance from the ionizing stars in cm. 
Ionization stages are identified by their color: dark blue (V), cyan (IV), green (111), orange (11), red (1). See Fig. 4 for the corresponding electron 
temperature and density profiles. 



underabundant in the neutral gas with respect to the ionized gas, 
although with a large uncertainty, with 5Hi(Si) ~ 0.71 + 0.42. 
Figures A. 2,A. 3 show that the line profiles are not consistent 
with the H ii region abundance. Silicon can be affected by de- 
pletion on dust in the neutral phase, but this turns out to be 
a negligible effect (Sect. 6.7). Furthermore, although a fraction 
of Si II could originate in the diffuse ionized gas (Sect. 4.2), it 
would only increase the discrepancy between the silicon abun- 
dance in the two phases. The Si/S ratio in the ionized gas is solar 
([Si/S]= -0.08+0.30), while it is somewhat sub-solar in the neu- 
tral gas ([Si/S]= -0.60 + 0.32). Assuming that [Si/S] should also 
be ~ in the neutral gas, we consider that our silicon abundance 
might be somewhat underestimated. This hypothesis seems to be 
confirmed in Fig. 11, which shows that our observed [Si/O] lies 
below the predicted value from massive star yields (Woosley & 
Weaver 1995). 

The oxygen abundance is widely used as a metallicity tracer 
However the O i lines in the COS spectra are either too weak 
to provide a useful upper limit or too strong to provide a re- 
liable column density determination (Sect. 3.3). To circumvent 



this issue, we can use silicon or sulfur as a replacement. Sulfur 
is a better candidate because the S ii lines are not prone to sat- 
uration. Bowen et al. (2005) find that S ii traces well H i in a 
diffuse slab over a wide range of ionization parameter (see also 
Sect. 6.1). Sulfur is thus a reliable substitute for oxygen and a 
useful metallicity tracer Assuming that the S/O ratio in the ion- 
ized phase of IZw 18 (-1.64) also holds in the neutral phase, 
we find (Os/H)a! 7.03 ± 0.15. This value is barely in agree- 
ment with the O/H ratio measured in the ionized gas, leading 
to (5hi(Os) ~ 0.18 + 0.16. Our indirect determination of the oxy- 
gen abundance lies between the FUSE value of L04, who found 
similar O/H between the two phases (although with large error 
bars), and the value of A03, who, despite using the same spec- 
trum as L04, found 6m(0) ~ 0.54 (Sect. 3.3). 

6.3. The odd-Z elements 

Although Al and P are not c-elements, they are thought to be 
produced the same massive stars that make oxygen. No signifi- 
cant amount should be produced in the explosive phases and the 
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Fig. 10. The abundance of each element is plotted in the H i region 
and in the H ii region. Elements are ordered by increasing condensation 
temperature (temperature for which 50% of an element is removed from 
the gas phase). The horizontal lines show the metallicity for each phase 
based on [O/H]. For the H i region, S is used to estimate O (Sect. 6.2). 



yields are expected to be metallicity dependent because of the 
odd-even effect (Woosley & Weaver 1995). Aluminum is pro- 
duced by C and Ne burning while phosphorus is produced dur- 
ing Ne burning. Lebouteiller et al. (2005) found that P n and O i 
column densities trace well each other over a wide range of envi- 
ronments as long as ionization corrections are not a major issue. 
Depletion on dust grains affects P more than O (Jenkins 2009), 
so part of the dispersion observed by Lebouteiller et al. (2005) 
could also be due to a differential depletion factor 

In this section, we discuss the abundances of P and Al rel- 
ative to ff-elements since there is no measurement of P or Al 
abundance in the ionized gas. Assuming a solar P/0 ratio we es- 
timate (Op/H)~ 6.86 + 0.33. This value agrees within errors with 
the oxygen abundance (Os/H)* 7.05 ±0.15 (Sect. 6.2). Figure 1 1 
shows that the expected yields from massive stars are also close 
to the solar P/O ratio. 

The Al/S ratio in the neutral gas is remarkably lower than the 
solar value, with [Al/S]= -0.74 ± 0.32. Aluminum abundance is 
determined using a single Al ii line that is unfortunately satu- 
rated. However, the saturation of the Al n line is not as strong as 
C n or Si ii lines, so we consider that saturation cannot explain 
the entire difference. It is also unlikely that there is a significant 
contribution from Al i given its extremely low ionization poten- 
tial of 5.99 eV. The low [Al/S] ratio could be due to a relatively 
stronger depletion of Al on dust grains as compared to S. The 
depletion of Al in the Local Interstellar Cloud is ~ 0.9 dex rela- 
tive to O in the WNM (Kimura et al. 2003). Such a factor would 
be more than enough to reconcile our Al/S ratio with the solar 
ratio. However, we investigate further the hypothesis of dust de- 
pletion in Sect. 6.7 and find that depletion alone cannot explain 
the low [Al/S] ratio. The Al/S ratio in the neutral gas, although 
sub-solar, is in fact consistent with expected yields from massive 
stars (Fig. 11). We therefore conclude that aluminum in the H i 
region is not deficient with respect to the a-elements. 

6.4. Carbon 

Carbon is produced in intermediate- and high-mass stars by He 
burning. During the early stages of a galaxy evolution, carbon 
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Fig. 11. Abundance ratios calculated from the stellar yields of Woosley 
& Weaver (1995). Yields are plotted for two mass ranges, 1 1 - 25 Mq 
(open symbols) and 30 - 40 Mq (filled symbols), and two metallici- 
ties, 0.1 Zq (diamonds) and 0.01 Zq (squares). The grey stripe shows 
the observed values. The horizontal black line shows the average ratio 
observed in the H ii region of BCDs (Izotov et al. 1999). 



is released mostly by the same young stars that make oxygen, 
so it is expected that C/O traces the yield ratio from massive 
stars. As galaxies evolve, older stellar population will contribute 
to the total carbon abundance, and a rise in C/O is expected (e.g., 
Izotov & Thuan 1999). 

In the ionized gas of IZw 18, [C/0]= -0.40 ±0.17, which is 
somewhat lower than more metal-rich BCDs, suggesting that C 
mostly comes from massive stars, as proposed in Izotov & Thuan 
(1999). Note that the C/O measured by Izotov et al. (1999), 
-0.77 ± 0.10, is slightly lower than both the value in Garnett 
et al. (1997), -0.63 ± 0.10, and than the value in Pequignot 
(2008), ~ -0.66. Garnett et al. (1997) found that C/O in IZw 18 
is larger than expected from massive star nucleosynthesis and at- 
tribute the discrepancy to C enrichment by a lower mass stellar 
population. 



18 



Lebouteiller et al.: Chemical enrichment and physical conditions in I Zw 18 



The carbon abundance in the H i region is somewhat lower 
than in the H n region. Using sulfur to trace oxygen in the neutral 
gas (Sect. 6.2), we find [C/Os]= -0.77 ± 0.33, which falls some- 
what below, although compatible with, the ionized gas value. We 
also note that [C/Os] appears to be compatible with expected 
massive stars yields (Fig. 11). 

6.5. Nitrogen 

The nucleosynthesis of nitrogen is not as well understood as for 
carbon or oxygen. Nitrogen can be produced in the CNO cycle 
by any stars more massive than 1.5 M©. The observation of a 
N/O plateau at low-metallicity in BCDs (e.g., Izotov & Thuan 
1999; Nava et al. 2006) leads to the conclusion that nitrogen is 
produced as a primary element, i.e., using the C and O produced 
in situ (since little C is available when the stars form). At higher 
metallicities (12 + log(0/H) > 8.3), secondary nitrogen pro- 
duction (from C and O from older generations) dominates. The 
origin of the N/O scatter at low-metallicity is still debated. If 
it is dominated by measurement errors, then nitrogen, like oxy- 
gen, is likely produced in short-lived massive stars (Nava et al. 
2006). The spread could also be intrinsic, due to the delayed 
injection of N from intermediate-mass stars. The observed scat- 
ter is theoretically explained either by effective yield variations 
and/or different star-formation history (Henry et al. 2000, 2006; 
Meynet & Maeder 2002b). The results of Henry et al. also show 
that nitrogen production in massive stars is possible. 

The N/O ratio in the ionized gas of I Zw 18, -1.61 ± 0.06, 
is compatible with the low-metallicity BCD plateau located at 
« -1.60 (Izotov & Thuan 1999). In the neutral gas, nitrogen is 
somewhat deficient, with 5hi(N) » 0.39 ±0.11. This translates 
into a N/O ratio of N/Os = -1.84 ± 0.19 which is only some- 
what lower than the ionized gas ratio. Following the N/O plateau 
paradigm, we conclude that the ISM was enriched on large scales 
by the same stellar population. Assuming the lower N/O in the 
H I region is genuine, it could be that intermediate-mass stars 
have contributed to the nitrogen enrichment in the H i region, or 
that nitrogen abundance was enhanced by massive stars in the 
H n region (see also Sect. 7.2). 

6.6. The iron group elements 

The iron group elements are particularly important because they 
could be released in the ISM on longer timescales than a- 
elements. Iron is produced partly in massive stars and in SNe la. 
A useful quantity to disentangle the two origins is [a/Fe]. The 
[a/Fe] ratio should increase with the progenitor mass (Woosley 
& Weaver 1995). If star formation is stopped, SNe 11 contribute 
no more to the Fe enrichment and SNe la could become impor- 
tant. Our Fe abundance is the same as the value determined by 
L04 using FUSE, and both values hint at a lower Fe abundance 
in the neutral gas as compared to the ionized gas. However, our 
value remains compatible with the ionized gas value within er- 
rors, with 5Hi(Fe) « 0.24 ± 0.28. We find that [Os/Fe] in the 
neutral gas is 0.45 ± 0.25, as compared to 0.37 ± 0.21 in the 
ionized gas. It is thus likely that the enrichment from SNe la, 
if any, is identical in the H i phase and in the H n phase. This 
value of [Os/Fe] lies between (solar O/Fe ratio) and the value 
in Galactic halo stars (Barbuy 1988) and in metal-poor stars in 
the solar neighborhood (Jonsell et al. 2005), suggesting an early 
enrichment by high-mass stars (see Izotov & Thuan 1999). 

Nickel and manganese both belong to the iron group. The de- 
pletion strength of Ni is similar to that of Fe while Mn is less de- 



pleted than Fe by a factor < 0.5 dex (Jenkins 2009). Manganese 
is mostly made in explosive Si burning while nickel is pro- 
duced mostly by neutron capture on the iron-peak elements in 
SNe la. It should be noticed that Ni and Fe are among the most 
heavily depleted elements in the ISM (e.g.. Savage & Sembach 
1996). We discuss the influence of depletion on dust grains in 
Sect. 6.7. The Ni/Fe ratio in the neutral gas of IZw 18 is solar, 
with [Ni/Fe]= -0.05 ± 0.23. We note that |Ni/Fe] remains essen- 
tially null over a wide range of metallicities in nearby stars (also 
in DLAs; Izotov et al. 2001b) while at the same time [a/Fe] in- 
creases at low-metallicity (Jonsell et al. 2005). This implies that 
Ni and Fe evolve in paraUel, and it is therefore not surprising 
that we also find a solar ratio in IZw 18. Since we already sug- 
gested that iron might be produced mainly in massive stars, the 
same result should apply to nickel. Concerning manganese, its 
abundance is uncertain, and we measure a Mn/Fe ratio slightly 
above the solar ratio, with [Mn/Fe] = 0.32 ±0.51. 

6.7. Depletion on dust grains 

We observe a small abundance discontinuity between the H i and 
the H II region (Sect. 6). Could part of the discontinuity be ex- 
plained by a differential depletion factor on dust grains? Volatile 
elements such as sulfur and chlorine provide us with a way of 
estimating the fraction of the underabundance in the neutral gas 
that is not due to depletion on dust. Unfortunately, CI is expected 
both as CI I and CI ii in the neutral phase so our upper limit 
on CI abundance is not useful. On the other hand, sulfur is our 
best metalUcity estimator since little or no S is expected to be 
locked in dust. We found that S/H is at least ~ 15% lower than 
in the ionized gas, implying that depletion alone cannot explain 
the lower metallicity of the neutral gas. 

We then examine whether depletion could explain the abun- 
dance discontinuity observed for the refractory elements. For 
this, we calculate the mass of each element in the H i region 
locked in dust that would be required to explain the observed 
abundance discontinuity between the H i and the H n region. 
Given the small mass fraction of ionized gas in IZw 18 (e.g., 
Petrosian et al. 1997), we expect most of the dust mass to lie in 
the H I region and we therefore ignore the depletion in the H n 
region. Dust extinction toward IZw 18 is mostly Galactic, and 
the dust mass in IZw 18 is exceptionally low. Using data from 
the Spitzer Space Telescope, GaUiano et al. (2008) found an up- 
per limit of < 2000 M0, while Herschel data constrains for the 
first time the dust mass, with 300;^™[5 M© (Remy et al. in prep.). 
We assume that the metals are well mixed in the H i region and 
express the mass locked in dust mi{X) as: 

md(X) = (/-!) mg(X) ~ (/ - 1) A(X) (X/H)obs Mm, d D 

where m^iX) is the mass locked into dust, OTg(X) is the observed 
abundance in the neutral gas, / is the total mass (dust-i-gas) di- 
vided by the mass in the gas, A(X) is the atomic mass, (X/H)obs 
is the observed abundance ratio, and Mm is the H i mass. The 
total H I mass in IZw 18 is Mri = 6.9 x 10^ Mq, 38% of which 
(2.6 X 10^ Mq) coming from the main H i component associated 
with the star-forming regions (Lequeux & Viallefond 1980; van 
Zee et al. 1998). 

The mass of each element in the gas phase is tabulated in 
Table 9. We also fist the mass required in the dust phase to reach 
the corresponding abundance in the H 11 region. For Al, P, and 
CI, whose abundances are not constrained in the ionized gas, we 
estimate the mass required to obtain a solar proportion with oxy- 
gen, i.e., [X/0]= 0. The mass of elements locked into the dust 
can be then compared to the dust mass measured by Herschel. 
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Table 9. Mass estimates of elements locked in dust. 



Element mg(Xy 


Am* 


(Mo) 


(Mo) 



N 

Os 

Al 

Si 

P 

S 

CI 

Fe 



■^^-^-160 

59+" 

4g70+1920 
^"'"-1365 

6 7+* '' 

" -3 3 

<4.7 

074+220 



59-86 
- 2100 
46-56 
210 - 450 
0.4 - 3.4 
33-95 
<4.4 
0-275 



Mass of element in the gas phase of the H i region (observed). 
Mass required in dust to reach the H n region abundance, i.e., 
(5hi(X) = 0. For Al, P, and CI, for which there is no H n region 
reference, we provide the mass required to reach the metallicity 
calculated from O, [X/H]= -1.64. 



Since we ignore the dust in the H 11 region, the dust mass de- 
rived from Herschel is an upper limit to the available dust in the 
H I region. We find that depletion alone cannot explain the un- 
derabundance of C and Si in the neutral gas. Even though the 
required mass of N, Al, and S into dust is below the observed 
dust mass, altogether they represent at least half of the available 
dust, requiring unlikely dust compounds composition. We thus 
conclude that depletion on dust cannot explain the observed un- 
derabundances in the neutral gas. 

Since there is little dust available to remove elements from 
the gas phase, we also find that the abundance ratios are not 
aflfected by depletion on dust. As an illustration, there is not 
enough dust to reconcile O/Fe with the solar ratio (Sect. 6.6). 
Similarly, we find that the low Al/S and Si/S ratio (Sects. 6.3, 
6.2) cannot be explained by aluminum and silicon locked in dust. 

7. Discussion 

The H I region of I Zw 18 is not pristine; it has already been en- 
riched in metals. We also find that all the elemental abundances, 
in particular for Al, S, P, Fe, and Ni, are consistent with produc- 
tion in massive stars (Sect. 6). Our best tracer of the enrichment 
in the H i region by massive stars is sulfur (Sect. 6.2). By using 
the solar sulfur abundance from Asplund et al. (2009), we es- 
timate the metalUcity of the neutral phase to be « l/46Zo, as 
compared to s; 1/317© for the H 11 region. Depletion on dust 
grains should have a negligible impact in both phases consider- 
ing the low dust mass in the galaxy (Sect. 6.7). We now discuss 
the enrichment of the H i and H 11 regions. 

7.1. Abundance discontinuity H i-H n 

The ionized gas in IZwl8 is more metal-rich than the neu- 
tral gas, with abundances higher by a; 0.18 - 0.77 dex depend- 
ing on the element. Our results are a priori compatible with ei- 
ther (1) the H II region being enriched by the present or recent 
star-formation episode(s) or (2) the presence of extremely low- 
metallicity clouds in the H i region diluting the observed abun- 
dances. 

The elements synthesized and released by the current burst 
are expected to lie at first in a hot phase (10* K) and mixing 
timescales could be as long as 10^ yr. Hot gas outflows have been 
observed through the O vi FUV doublet in starburst galaxies, in- 
cluding a few dwarf galaxies (e.g.. Grimes et al. 2009). However, 
there is no evidence that I Zw 18 developed such an outflow, even 



though the non-detection of the O vi doublet could be due to 
the low oxygen abundance in this object (Heckman et al. 2002). 
Furthermore, Chandra X-ray observations of SBS 0335-052W, 
SBS 0335-052E, and IZw 18 also failed to reveal any sign of hot 
gas breaking out from the stellar body implying that the galactic 
fountain scenario may not be valid in BCDs (Thuan et al. 2004; 
see also Sect. 7.3). 

As first hypothesized by Kunth & Sargent (1986), the H 11 
region could be self-enriched, with a fraction of the newly pro- 
duced heavy elements dispersing on small spatial scales and 
mixing on short timescales (few Myrs). This would naturally ex- 
plain the H i-H n region abundance discontinuity. Legrand et al. 
(2000) find that abundances are remarkably uniform in the ion- 
ized gas of IZw 18, with a dispersion of < 0.1 dex on scales of 
several hundreds of parsecs, which seems to go against the fast 
and local enrichment hypothesis in the H 11 regions. We recall 
that our H i region abundances are measured toward the NW re- 
gion, and compared to the H 11 region abundances in that region. 
The NW region is where the youngest stars are located, which in 
principle provides us with the best chance of observing a local 
enrichment from the recent star-formation episode(s). The most 
significant abundance discontinuities H i-H 11 we calculate are 
for N and Si, with 5hi(N)= 0.39 + 0.14 and 5Hi(Si)= 0.71 + 0.42 
respectively (Sects. 6.5, 6.2), which is significantly larger than 
the abundance dispersion across the galaxy in its ionized gas. 
Although nitrogen abundance might be enhanced by very mas- 
sive, rotating, stars in the H 11 region (e.g., Meynet & Maeder 
2002a; Hunter et al. 2008; Maeder & Meynet 2012), local en- 
richment by massive stars does not explain the abundance en- 
hancement of the other elements in the H n region. 

Another possibility to explain a metal deficiency in the neu- 
tral gas is to invoke the presence of extremely metal-poor regions 
along the lines of sight, such as the low-metallicity high- velocity 
clouds infalling on the Milky Way. The H i envelope in IZw 18 
extends as much as ~ 13kpc away from the star-forming re- 
gions (Lelli et al. 2012), and it is possible that the - presumably 
lower density - outskirts are not as enriched in heavy elements 
as the immediate surroundings of the H n regions. In order to 
estimate the mass of this relatively metal-poor H i gas, we ten- 
tatively assume that the H i medium can be separated into two 
components, one (hereafter Hl-sf) is spatially associated with 
the star-forming region and another component (Hl-ex) is more 
extended and possibly less enriched. In the following, we neglect 
the size and mass of the H n region. We further assume that the 
mixing of metals is uniform in the vicinity of the H 11 region, i.e., 
abundances are similar in the H 11 region and in Hl-sf. The metal 
abundance in the H i region can then be written as 



1 ^" /x fv\ px 

(X/H)hi = f (X/H),f 



(12) 



1+/h 

And the underabundance (Jhi(X) in the neutral gas becomes: 

5hi(X) - - log [( 1 - /h) + /x /H,exl , ( 1 3) 

where fx - (X/H)ex/(X/H)sf is the relative abundance ratio of 
a given element between Hl-ex and Hl-sf and /h is the fraction 
of the total H i along the line of sight in the Hl-ex component. 
Figure 12 shows the variations of (5hi(X) with /h and /x. 

Taking sulfur as a metallicity tracer (X -S), it can be seen 
that if Hl-ex is pristine (/s = 0), the observed sulfur deficiency 
is compatible with « 10 - 45% of hydrogen into Hl-ex along 
the line of sight. Assuming a spherical distribution with constant 
hydrogen volume density, this represents a factor [1 - (1 - /h)^] 
in mass, translating to ~ 25 - 80% of the total H i mass in Hl- 
ex. In other words, about half the mass of the H i envelope in 
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Fig. 12. The metal deficiency in the H i region dm is plotted as a func- 
tion of the fraction of hydrogen present in the Hl-ex component defined 
as the extended H i component that may be less enriched in metals as 
compared to Hl-sf (the H i component associated to the star-forming 
regions). We assume that the Hl-sf component has the same metallicity 
as the H n region and calculate 6m for several values of the sulfur abun- 
dance ratio between Hl-ex and Hl-sf. The horizontal dashed line and the 
grey stripe indicates illustrate the sulfur deficiency, ^hi(S) = 0.18±0.1 1. 



IZw 18 could be pristine. This scenario is therefore not compat- 
ible with a small number of low column density H i clouds such 
as the high-velocity clouds observed in the Milky Way. If the 
observed H i-H ii region abundance discontinuity is caused by 
the presence of pristine gas, the latter should represent a signifi- 
cant mass fraction of the galaxy (whose baryonic mass budget is 
dominated by H i gas). 

We can then assume that the Hl-ex component has also been 
somewhat enriched (/s > 0). Considering that the H i mass in 
Hl-ex is diluted over a larger volume, the aperture fraction is 
probably very small, and we can reasonably assume that the frac- 
tion of hydrogen in the Hl-ex component is less than 75%. For 
such a value of /h, the observed sulfur deficiency is compati- 
ble with the sulfur abundance in Hl-ex being =s 40 - 80% of the 
value in Hl-sf, and thus also of the value in the H ii region. This 
would imply a metallicity as low as ~ 1/702© in the outskirts 
of the galaxy. These results could put galaxy evolution scenar- 
ios to the test, with in particular the hypothesis of a minimal 
enrichment level at a given redshift (Sect. 7.3). Does a metallic- 
ity of ~ 1 /70 Zq correspond to the minimal enrichment at red- 
shift z = 0? Deep H i surveys such as ALFALFA (Giovanelli 
et al. 2005), and later SKA, should identify large numbers of 
optically-faint H i clouds in the Local Universe that will test the 
hypothesis of a minimal enrichment level at z = 0. We also wish 
to emphasize that the hypothetical presence of metal-free gas has 
important consequences on the interpretation of the integrated 
D/G ratio, in particular if the latter does not scale linearly with 
metallicity. 

7.2. Comparison with the sample of BCDs 

We now compare our results on IZw 18 to the sample of BCDs 
whose neutral gas has been investigated with the FUSE tele- 
scope. First, we examine the N/O abundance ratio as it poten- 
tially holds clues as to the stellar population responsible for en- 
riching the ISM. Nitrogen can be synthesized in a wider range 
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Fig. 13. Top - The N/O ratio is compared between the H i and H n 
region for the sample of BCDs studied with FUSE, and for IZw 18 
with COS. The labels are defined as la: IZw 18 (L04), lb: IZw 18 
(A03), Ic: IZw 18 (this study, with S as a tracer of O; Sect. 6.2), 2: 
NGC 1705 (Heckman et al. 2001), 3: SBS 0335-052 (Thuan et al. 2005), 
4: Mark 59 (Thuan et al. 2002), 5: NGC 625 (Cannon et al. 2005), 6: 
IZw 36 (Lebouteiller et al. 2004), 7: Pox 36 (Lebouteiller et al. 2009), 
8: NGC 604 (Lebouteiller et al. 2006), 9: SBS 1543+593 (Bowen et al. 
2005; Schulte Ladbeck et al. 2004). Bottom - Ratio of the average metal 
abundance [M/H] =< [X/H] > between the H i and H ii regions. [M/H] 
is calculated two ways, averaging all elements observed in each phase 
(small diamonds), and averaging only the elements observed in both 
phases (large diamonds). 



of stellar masses than oxygen, and rapidly rotating massive stars 
can potentially enrich the ISM on short timescales (Sect. 6.5). In 
contrast, oxygen and the a-elements are produced exclusively in 
massive stars and are released within a few Myrs through SNe 
(Sect. 6.2). Figure 13a compares the N/O ratio in the H i and in 
the H II region of BCDs. The N/O ratio is similar in both phases 
for the highest metallicity objects. There might be a hint of a 
lower N/O in the H i region for the lowest-metallicity objects 
(I Zw 18, SBS 0335-052, and IZw 36), but the result is inconclu- 
sive within error bars. We conclude that the stars producing N 
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have the same mass range in both phases, which argues for an 
enrichment of the ISM before the present starburst episode. 

The hypothesis of relatively metal-poor outskirts (Sect. 7.1) 
can also be put to the test in the BCD sample. Lebouteiller et al. 
(2009) summarized the FUSE results and showed that the metal- 
licity in the H I region in BCDs was equal (IZw 18, SBS 0335- 
052) or lower (I Zw 36, Mark 59, Pox 36, NGC 625, NGC 1705) 
than that in the H n region. A trend might exist, with the discon- 
tinuity being essentially null for the lowest metallicity objects, 
and then constant around ~ Idex for metallicities > 1/6 Zg 
(value in the ionized gas). However, if nitrogen is used as a 
metallicity tracer instead of oxygen, the results could also be 
interpreted as a constant discontinuity independent on metallic- 
ity. The lowest-metallicity object, IZw 18, is key to conclude 
on the ubiquity of the abundance discontinuity H i-H n region. 
While the result for IZw 18 was at the time inconclusive be- 
cause of the discrepancy between A03 and L04, our study sug- 
gests that the abundance discontinuity in I Zw 1 8 is rather small, 
smaller than what is observed in other BCDs. We now revisit 
the findings of Lebouteiller et al. (2009) by including this new 
measurement, and by calculating the average metal abundance 
[M/H] [X/H] > in the H i and H ii regions (Fig. 13b). The 
choice of averaging the abundances is motivated by (1) possible 
systematic problems in the abundance determination of some el- 
ements (hidden saturation, stellar contamination, atomic data...), 
and (2) the fact that there is so far no clear evidence of a different 
nucleosynthesis origin of elements in both phases. Figure 13b 
shows that the abundance discontinuity is around 0.7 dex for 
most objects, with no clear correlation as a function of metal- 
licity. This imphes that the chemical enrichment H i and H n 
region is disconnected in most objects, with the H i region sys- 
tematically more metal-poor than the H n region. Applying the 
results of Sect. 7.1 to the sample of BCDs, Figure 12 shows that 
for ^Hi ~ 0.7 dex, the Hl-ex component (extended H i enve- 
lope) would have a metallicity between 0% and 30% that of 
the Hl-sf component (associated to the H ii region), and would 
represent more than 75% of the H i mass. This result is rele- 
vant in particular for SBS 0335-052, which, despite a metallicity 
close to I Zw 18, shows a discontinuity of « 0.7 dex (Figure 13b). 
For SBS 0335-052, the H i envelope could be as metal-poor as 
~ 1/150Zq. New abundance determinations in this object by 
COS should bring a definitive answer (James et al. 2013). 

7.3. Enrichment of the H i region 
7.3.1. Cosmic evolution 

Despite the weak discontinuity between the abundances of the 
H I and H n regions of IZw 18, the ISM is fairly well mixed on 
kiloparsec scales, with the metallicity of the neutral gas at least 
50% of that in the ionized gas (most hkely value of 66%). 
How was the ISM, and in particular the H i region, enriched? 
Considering the relatively uniform abundances, at least half of 
the metallicity in the H i region was presumably set before the 
current starburst episode. A major uncertainty concerns the gen- 
uine age of the galaxy as compared to its chemical age. The latter 
is probed by the metallicity parameter, which is determined by 
a complex interplay between successive star-formation episodes 
and external gas flow exchanges (outflow of enriched material 
and infall of metal-poor gas). 

Although the age of the oldest stars in IZw 18 is still sub- 
ject to debate (see Papaderos & Ostlin 2012), stars as old as 
the Hubble time (Aloisi et al. 2007; Contreras Ramos et al. 
201 1) might be present, suggesting that 1 Zw 18 is not genuinely 



young, and therefore that its present-day metallicity is the re- 
sult of the mechanisms mentioned above over cosmic times. For 
this reason, the comparison between IZw 18 and high-z objects 
is relevant and instructive. The intergalactic medium (IGM) was 
already enriched to metallicities similar to that of the H i re- 
gion of IZw 18 at z ~ 1.6-2.9 (Telfer et al. 2002). Hence 
it is possible that IZw 18 was assembled from a collection of 
low column-density clouds in the IGM a few Gyrs after the Big 
Bang, and experienced little enrichment since then. The metal- 
licity of the H I envelope of I Zw 1 8 is also typical of DLAs at 
z > 2 (Rafelski et al. 2012), again compatible with a slow chem- 
ical enrichment. The lack of strong constraints on the mass of 
DLAs makes the comparison with I Zw 1 8 relatively uncertain 
within the chemical downsizing scenario in which galaxies with 
different mass have different star-formation efficiencies (e.g.. 
Brooks et al. 2007). Nevertheless, the lack of nearby galaxies 
as metal-poor as high-z DLAs (~ 1/100 Z©; e.g., Prochaska & 
Wolfe 1999; Rafelski et al. 2012) implies that nearby BCDs, in- 
cluding IZw 18, have experienced a slow, minimal, enrichment 
over several Gyrs. Legrand et al. (2000) suggested that IZw 18 
went through a mild continuous star-formation rate (SFR), ~ 
10""* Mq yr"', over several Gyrs as opposed to a series of violent 
short bursts separated by relatively longer quiescent eras. Metal 
enrichment in BCDs through a series of violent star-formation 
episodes is an unhkely process as galaxies with difl'erent star- 
formation histories would result in diff'erent metalUcities, thus 
including a class of nearby galaxies with no starburst episodes 
and an extremely low-metallicity, which is not confirmed by ob- 
servations. 

The relation between stellar mass and metalUcity in galaxies 
(e.g., Tremonti et al. 2004) suggests that star formation proceeds 
slowly at low metallicity. However, the mass-metallicity rela- 
tion can also be understood as a lower star-formation efficiency 
at low mass, with a relatively more efiicient ejection of enriched 
material, or by a metallicity increase through star-formation trig- 
gered by infalling gas. Mannucci et al. (2010) showed that a 
more fundamental relation exists between stellar mass, metal- 
licity, and SFR, where downsizing, outflows, and infalls are the 
major actors. 

Accretion of inflowing cold gas triggers a transient phase in 
which (1) metallicity is diluted by the injection of significant 
amounts of low-metalhcity gas, (2) the SFR increases through 
the Schmidt-Kermicutt law, and (3) metaUicity increases again 
as star formation proceeds. Although H i companions are often 
observed around star-forming dwarf galaxies (e.g., Taylor 1997; 
Ramya et al. 2009; Pustilnik et al. 2001), their accretion and the 
consequence on star formation is not definitive (e.g., Telles & 
Maddox 2000). Gas inflow is expected to play a more important 
role at high-z. The cosmological inflow from the cosmic web 
could have been so high at high-z that star-formation could not 
catch up in the least massive haloes like IZw 18 (Krumholz & 
Dekel 2012). This metallicity-dependent quenching would result 
in a growth of the H i envelope, the latter still being observed 
today. Only at z < 2 could a low-mass galaxy like IZw 18 start 
to build up its metallicity. Such a scenario remains compatible 
with the presence of old stars in I Zw 18. 

Alternatively (or possibly in parallel), galactic winds could 
expel most of the metals into the IGM through multiple super- 
novae explosions, effectively keeping the ISM metallicity to low 
levels (e.g., Recchi et al. 2004). The most active star-forming 
regions would be responsible for those winds. A major uncer- 
tainty remains however on the fate of the metals observed in 
galactic outflows, since the superbubble expansion might de- 
celerate, and the outflowing material might cool down and re- 



22 



Lebouteiller et al.: Chemical enrichment and physical conditions in I Zw 18 



turn to the ISM instead of escaping the gravitational well (e.g., 
Tenorio-Tagle et al. 1999; Mac Low & Ferrara 1999; Rieschick 
2003). Despite the current starburst episode in IZw 18, there is 
no evidence of outflowing hot gas (Sect. 7.1). Furthermore, we 
note that outflows would play a more important role in galaxies 
whose star-formation history is dominated by a series of violent 
bursts, which has probably not been the case in I Zw 1 8 (Legrand 
et al. 2000). 

In summary, although gas flows play a major role in the cos- 
mic evolution of galaxies, their efi'ects in I Zw 18 are not yet fully 
constrained and understood. A combination of metal-poor gas 
infall in the early stages of the galaxy, and a chemical downsiz- 
ing could be responsible for the extremely low-metallicity ob- 
served today in the H i and H n regions. Finally, our results do not 
allow putting constraints on the reasons of the current starburst 
episode in IZw 18, which has been ~ 10-100 times relatively 
stronger over the past ~ 50 Myr (Aloisi et al. 1999; Legrand et al. 
2000). The sudden onset of star formation could have been trig- 
gered by gas compression from merging events. The disturbed 
H I morphology suggests that several small structures might have 
been interacting (Lelli et al. 2012). Such structures could have 
been gas-rich dwarf galaxies or sub-DLAs. From the radial ve- 
locity of the neutral metals observed with COS, we did not find 
any evidence of infalUng material on the NW region itself; the 
H I gas is, if anything, outflowing from the H n region with a 
velocity of ~ 20kms"' (see also A03; L04). However, Lelli 
et al. (2012) found a radial inflow/outflow motion on the order of 
~ 15 km s"' on larger scales and concluded that interactions or 
mergers are the most Ukely explanation for the current starburst 
episode. In general, the FUSE results for other BCDs show blue 
shifts of at most 10 - 20kms"' (e.g., Lebouteiller et al. 2009; 
Cannon et al. 2005), confirming that slow radial flows of H i gas 
might be frequent. 

7.3.2. Local conditions for star formation 

We now discuss the star-formation process itself, and the influ- 
ence of metallicity in the ability of clouds to form stars in an 
environment as metal-poor as IZw 18. Our present study puts 
strong constraints on the metallicity of the H i region, which 
constitutes the bulk of the mass of the galaxy. The H i region is 
thus key to understand the past star-formation efficiency. 

First, the fact that stars of all masses have formed for several 
Gyrs has some implications for the early metallicity evolution of 
the galaxy. Brormn & Loeb (2003) showed that when the metal- 
licity reaches more than 10"* - 10"^ Zg, the cooling by the far-IR 
fine-structure lines [C n] 157 /urn and [O i] 63 /im enables cloud 
fragmentation and therefore low-mass star-formation. No dust 
or molecules were considered for this calculation and we note 
that H2 could provide the most efficient cooling at even lower 
metallicity (e.g., Omukai et al. 2005). However, the lack of dif- 
fuse H2 in IZw 18 (Vidal-Madjar et al. 2000) suggests that H2 
cooling is not significant in our practical study, except if present 
in small dense clumps where H2 might be more abundant (see 
Thuan et al. 2005). We thus get the picture of IZw 18 forming 
stars of all masses only a few Gyrs after the Big Bang, when 
its metallicity reached > lO^^^Z©. The oldest stars in IZw 18 
could be the remnants of the star-formation episode following 
this early enrichment. The existence of a metallicity threshold 
for low-mass star-formation is however still debated and we note 
that the presence of dust to shield UV photons, and the D/G ratio, 
ultimately related to metallicity via the evolution of the dust-to- 
metal ratio, may be the most important parameters governing the 



star formation process (Schneider et al. 2012; Glover & Clark 
2012). 

The importance of the ISM metallicity is fundamental for the 
ability of clouds to form stars, and it may well regulate the SFR 
as long as external gas flows are unimportant. Although the ob- 
served SFR in BCDs is ~ 10"^ - lOMoyr \ this corresponds 
to the current starburst episode, and a much milder SFR could 
have dominated beforehand. For example, in IZw 18 the past 
SFR could have been as low as ~ lO'^'Mgyr"'; e.g., Legrand 
et al. 2000). It seems at first contradictory with Mannucci et al. 
(2010) who finds that for a given stellar mass, the SFR is anti- 
correlated with metallicity. The fundamental metallicity rela- 
tion found by Mannucci et al. (2010) is, however, biased toward 
emission-line galaxies, and a "pre- starburst" IZw 18 would not 
have been selected by their survey. Furthermore, as discussed in 
the following, it is possible that the extremely low-metalhcity 
of this galaxy has a major impact on the star-formation process 
itself. 

The temperature and density of the ISM are essential param- 
eters for star formation. The latter can be inhibited by the inabil- 
ity of the gas to cool down enough for removing thermal pressure 
and for creating gravitational instabihties. The warm tempera- 
ture in a metal-poor ISM is due to the paucity of efficient cooUng 
mechanisms, such as molecules or the far-IR fine-structure fine 
[C n], which dominates the cooling in the diffuse atomic clouds. 
Krumholz (2012) showed that the thermal equilibrium timescale 
is more important than the timescale for conversion of atomic 
to molecular gas in a low-metalUcity gas. Star formation could 
proceed before the medium becomes molecular when the metal- 
licity reaches below a few percents solar. The metallicity thresh- 
old above which thermal equilibrium is reached before a few 
free-fall timescales (once the gas becomes gravitationally unsta- 
ble) thus depends mostly on the abiUty of the gas to cool down. 
The H I region of IZw 18 is an ideal testbed for star-formation 
timescales in a low-metallicity atomic gas, since molecular gas 
has yet to be detected. The model in Krumholz (2012) highly 
depends on the heating mechanisms and their metalhcity depen- 
dence, and the uncertain heating mechanisms in the H i region 
of galaxies such as IZw 18 (Sect. 5.4) add to the uncertainty 
on the threshold metallicity and thermal equilibrium timescale. 
Following the method in Krumholz (2012), we use our determi- 
nation of the heating rate in the H i region, and an initial gas den- 
sity and temperature of ~ 10 cm~^ and ~ 1000 K, respectively. 
We find that thermal equilibrium is reached before a few free- 
fall timescales only for a metallicity greater than ~ I/IOOZq. 
A lower threshold metallicity would be obtained if we consider 
a significant clumping factor. Hence we see that widespread star 
formation is possible in the H i region of 1 Zw 1 8, currently hold- 
ing a metallicity of » 1/46 Zq. If such physical conditions pre- 
vail in objects as metal-poor as < l/lOOZg, the thermal equi- 
librium timescales are rather long and the abiUty to form stars 
ultimately depends on the survival of the clouds. IZw 18 might 
have kept a low metallicity because its ISM remained warm on 
long timescales. A major uncertainty remains on the presence of 
small cold clumps (cold neutral medium phase) in which star 
formation would be more efficient. The mixing factor of the 
warm neutral medium and cold dense clumps in the FUV absorp- 
tion lines is difficult to quantify. Moreover the ISM pressure in 
I Zw 18 might be too low to maintain a multi-phase ISM, so that 
the neutral gas we characterized in this study might be the only 
thermally-stable phase. The temperature from this phase could 
thus be a crucial parameter to understand the star-formation pro- 
cess on large spatial scales and over cosmic times. 
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8. Conclusions 

We present a comprehensive investigation of the H i envelope of 
IZw 18. The H i region holds several fundamental clues to the 
chemical evolution of galaxies. The metalUcity of the H i region 
is set by the successive star-formation episodes over the galaxy 
history, and is modulated by outflows of supernovae products 
and hypothetical inflows of quasi-pristine gas from the cosmic 
web. The abundance of metals in the H i region is also a critical 
parameter to understand the present and the future star-formation 
history, as metals provide most of the cooling in the H i gas reser- 
voir to facilitate and sustain the star-formation process. Our ob- 
jectives are to characterize the H i region by (1) its physical con- 
ditions and (2) its chemical abundances. We compare the abun- 
dances to those derived in the H ii region (from optical emission- 
lines), and discuss the enrichment history of the galaxy. 

The H I envelope of IZw 18 was observed toward the NW 
young stellar cluster in the FUV with the COS instrument on 
Hubble. The massive stars provide the FUV continuum over 
which absorption-lines from neutral species are superimposed. 
We modeled the stellar absorption in order to remove its con- 
tribution and isolate the interstellar absorption. We derived the 
column density of H i and of heavy elements (C n, C n*, N i, O i, 
Al n. Si n, P n, S n, Mn n, Fe n, Ni n). Models were constructed to 
compute the ionization structure of the nebula. Our models pre- 
dict a large ionized gas volume in the H n region, suggesting that 
the H I cavity observed at 21 cm toward NW is due to ionization 
by the stellar cluster as opposed to blow-out by stellar winds and 
supernovae. We study potential emission components in the H n 
region that could work out to reduce the apparent optical depth 
of absorption lines. 

The physical conditions in the H i region are drawn from the 
study of the C n* absorption-line. 

- The population of the fine-structure level of C"^ toward NW 
is compatible with collisions with H" atoms, implying that 
C n* arises in the H i envelope. 

- We determine observationally that the electron fraction in the 
H I region is less than ~ 2%; our models predict ~ 0.1%. 

- From C n* we derive a cooling rate that is larger than the 
models. An extra heating source is required in the models to 
reconcile them with the observations, leading to an elevated 
electron temperature in the H i region of ~ 1000 K. 

- The cooUng rate inferred from C n* is used to calculate the 
star-formation rate. Our value agrees with other determina- 
tions in the literature, although on the high-end. The D/G ra- 
tio used for this calculation is a factor of ~ 2000 times lower 
than the MiUcy Way value, i.e., a factor of ~ 60 times lower 
than what is expected if we simply scale D/G with metaUic- 
ity. 

We then infer the dominant ionization stages in the H i region 
from the models and calculate the chemical abundances in the 
H I region. 

- We find that sulfur is less abundant in the H i region as com- 
pared to the H II region, by a (log) factor of 0.18 + 0.1 1. 

- Nitrogen, and possibly carbon, are also underabundant in the 
H I region. 

- The P/S ratio is solar while Al/S is consistent with yields 
from massive stars. The Fe abundance is the same in the 2 
phases within errors. The [or/Fe] ratio is consistent with pro- 
duction in massive stars. 

- Depletion on dust grains is unlikely to explain the observed 
abundance discontinuity between the H i and H n regions. 
This is due to the low dust mass available. 



Finally, we discuss possible scenarios concerning the abundance 
difi'erence in the gas phases, and concerning the enrichment of 
the H I region. 

- Local enrichment in the H ii regions could explain the en- 
hanced N abundance but fails to explain other results, espe- 
cially for S. 

- The presence of relatively lower metallicity outskirts in the 
H I region is discussed. As much as 50% of the H i mass 
could be pristine. 

- The metallicity of the H i region is not null, it is compati- 
ble with the IGM metalUcity at z = 1.6 - 2.9. IZw 18 could 
be considered as an old DLA with little recent star forma- 
tion. The presence of old stars confirm that IZw 18 is not 
a genuinely young galaxy. The inflow of quasi-pristine gas 
from the cosmic web for a few Gyrs after the Big Bang 
could explain the massive H i envelope and the quenched 
star-formation in the early stages of the galaxy evolution. 

- The current onset of star formation could be a due the merg- 
ing of dwarfs or sub-DLAs, as suggested by the disrupted H i 
morphology. 

The sample of BCDs whose neutral gas has been investi- 
gated by FUSE and COS shows that a systematic discontinuity 
exists between the H i and H ii region abundances. IZw 18 is the 
object with the lowest discontinuity. Another BCD, SBS 0335- 
052, despite having a metalUcity similar to IZw 18 does show a 
significant discontinuity which could translate in the H i region 
being as metal-poor as 1/150 Z©. 
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Appendix A: Line fitting 
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Fig.A.l. Line profile fitting for the 11523 (top), 12028 (middle), and 11579 datasets (bottom). For display purposes, only the strongest lines 
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Fig. A.2. See Fig. A. 1 for the plot description. The unlabelled line on the far right is Si ni /i 1206.5 arising in IZw 18. 
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Fig. A.3. See Fig. A. 1 for the plot description. 
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Fig. A.4. See Fig. A. 1 for the plot description. 
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Fig. A.5. See Fig. A. 1 for the plot description. 
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Fig. A.6. See Fig. A. 1 for the plot description. 
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Fig. A.7. See Fig. A. 1 for the plot description. 
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